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the  material  parameters  needed  to  calculate  the  TLC  can  be  obtained  from 
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formability  of  certain  high  strength  aerospace  sheet 


alloys. 


Constant  extension  rate  tests  vs.  strain  rate  and  temperature  and  load 
relaxation  tests  vs.  temperature  have  been  performed  on  2024-0  aluminum  sheet 
and  annealed  T1-6At-4V  sheet  and  round  bar  tensile  specimens.  Load-elongation 
and  load-time  data  were  recorded  digitally  for  direct  determination  of  char- 
acteristic derivative  plots  such  as  dzncr/de  (v)  vs.  a and  o vs.  e where  o,  e, 
and  e are  the  stress,  non-elastic  strain  and  strain  rate  respectively.  Thus 
the  strain  hardening  coefficient  and  strain  rate  sensitivity  are  obtained  as 
functions  of  stress,  strain  rate,  and  temperature.  The  plastic  anisotropy  of 
the  sheet  materials  has  also  been  Investigated. 


FLC's  can  be  calculated  from  constitutive  equation  parameters  determined 
directly  from  the  derivative  plots.  From  the  dependence  y (o)»  it  was  found 
that  the  Swift  equation  models  the  strain  hardening  behavior  exactly  at  high 
strains  (e  > 0.05  for  2024-0  At;  e>0.02  for  T1-6AA-4V).  At  25°C,  the  strain 
rate  dependence  c (e)  could  be  modelled  by  an  equation  suggested  by  Hart  for 
low  homologous  temperatures.  For  2024-0  At  at  25°C,  the  strain  rate  sensi- 
tivity was  small  enough  to  be  neglected. 


For  2024-0  aluminum  at  25°C,  the  FLC  calculated  from  Hill's  theory  of 
localized  necking  (strain-rate  independent)  was  in  good  agreement  with  an 
experimental  in-plane  forming  limit  diagram  (FLD).  However,  both  the  calcu- 
lated FLC  and  in-plane  FLD  were  significantly  lower  than  the  standard  out- 
of-plane  FLD  determined  from  punch-stretch  tests.  For  2024-0  aluminum  at 
163°C  (325°F)  and  Ti-6At-4V,  an  FLC  based  on  a strain-rate  dependent  model 
will  be  required  for  comparison  to  experimental  limit  strains. 
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I.  INTRODUCTION 


1 . 1 Sheet  Metal  I'm  inability 

The  formability  of  sheet  metal  is  dependent  upon  certain  plastic  pro- 
perties of  the  material  identified  as  formability  indices.  These  indices 
affect  formability  in  two  distinct  ways.  On  one  hand,  they  determine 
the  limiting  states  of  strain  which  a sheet  can  provide  in  a forming  pro- 
cess without  an  intervening  failure.  In  addition,  they  affect  the  way 
strain  is  distributed  in  a sheet  during  forming  and,  hence,  partially  de- 
termine the  maximum  strains  that  will  be  imposed.  Hence,  it  is  necessary 
t.o  identify  and  determine  these  formability  indices  for  input  to  both  form- 
ability  and  process  models. 

Several  material  properties  are  widely  accepted  as  important  sheet 
formability  indices.  These  include  the  strain  hardening  coefficient  and 
strain  rate  sensitivity.  It  is  less  widely  recognized  that  these  indices 
can  show  a significant  dependence  on  strain-rate,  temperature,  and  some 
measure  of  plastic  hardness  due  to  prestrain.  It  follows  that,  in  order 
to  include  these  parameters  into  analytical  formability  or  process  models, 
it  is  important  to  have  an  analytical  material  model  or  constitutive  equa- 
tion which  can  be  used  to  express  these  dependencies. 

Recently,  Thomas,  Gegel,  and  Teutonico  [1]  presented  a systematic, 
approach  to  assessing  sheet  metal  formability  based  upon  analytical  mate- 
rial constitutive  equations  and  calculated  forming  limit  curves.  The  forming 
limit  curve  (FLC),  which  gives  the  locus  of  principal  strains  beyond  which 
failure  is  expected,  is  based  on  a specific  failure  mode  and  corresponding 
set  of  material  properties.  The  constitutive  equation  is  used  to  input 
these  properties,  including  their  dependence  on  strain,  strain-rate,  and 
temperature,  into  the  FLC  analysis.  When  combined  with  numerical  sheet 
metal  process  simulation  models,  an  integrated  CAD/CAM  program  for  sheet 
metal  forming  can  be  developed. 

1 . 2 Formability  Indices 

It  is  usually  assumed  that  the  most  important  failure  criteria  for 
sheet  forming  are  related  to  plastic  instability  since  plastic  instability 
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usually  precedes  fracture.  In  this  case  the  relevant  sheet  metal  forma- 
billty .indices  can  be  taken  to  be 

(1)  The  strain  hardening  coefficient  r = (n  n a / S c)  • where 

o is  the  flow  stress,  e is  the  non-elastic  strain,  and  e is  the 
non-elastic  strain  rate, 

(2)  The  strain  rate  sensitivity  v = ( 8 i.  n o / sm  c) 

(3)  One  or  more  history  dependent  parameters  describing  plastic 
hardness  which  evolve  during  deformation  and  which  also  partially 
determine  y and  v,  and 

(4)  Parameters  needed  to  describe  plastic  anisotropy  of  the  candidate 
sheet  material . 

The  functional  dependences  of  the  indices  (1)  and  (2)  on  stress  o and  strain 
rate  e and  the  index  list  (3)  will  depend  on  the  form  of  the  constitutive 
equation. 

The  indices  (4)  describing  plastic  anisotropy  are  taken  to  be  the 
plastic  strain  ratios  rQ,  rgQ,  and  r.g  measured  on  tensile  specimens  cut 
from  the  sheet  at  0°,  90°,  and  45°  to  the  rolling  direction  respectively. 

In  formability  calculations  which  depend  on  normal  anisotropy  only  [1],  an 
effective  r = 1/4  (ro+r90+2r45^  1S  used-  In  calculations  which  account  for 
planar  anisotropy  [2],  the  yield  function  depends  explicitly  on  all  three 
measured  strain  ratios. 

1.3  Scope  of  the  Program 

The  objective  of  this  program  has  been  to  determine  and  experimentally 
verify  material  analytical  models  or  constitutive  relations  applicable  to 
the  formability  and  processing  of  high  strength  aluminum  and  titanium  alloy 
sheet  materials.  The  constitutive  relations  will  be  used  in  related  Air 
Force  programs  on  modeling  sheet  formability  and  sheet  forming  processes. 

The  candidate  materials  selected  for  the  test  program  were  2024-0 
aluminum  and  the  titanium  alloy  Ti-6A£-4V.  The  2024  aluminum  Is  usually 
formed  In  the  "0"  or  soft  temper  and  then  heat  treated  to  high  strength. 
T1-6A&-4V  is  the  most  widely  used  titanium  alloy  for  aircraft  structural 
parts. 


2 


The  experimental  test  program  has  been  restricted  to  uniaxial  tensile 
tests.  These  Include  both  constant  extension  rate  tests  used  to  determine 
strain  hardening  behavior  and  load  relaxation  tests  used  to  determine  strain 
rate  sensitivity.  The  tests  have  been  conducted  on  commercial  sheet  mate- 
rials at  temperatures  representative  of  typical  forming  operations.  The 
plastic  anisotropy  of  the  sheet  has  been  considered.  A digital  data  acqui- 
sition system  has  been  developed.  Data  in  digital  form  Is  required  for 
accurate  determination  of  required  derivatives  and  for  convenient  mathe- 
matical modeling. 

The  uniaxial  tensile  test  configuration  has  been  selected  because  it 
provides  a uniform  stress  state  and  is  widely  used.  The  advantages  of 
being  able  to  characterize  the  formabllity  of  a sheet  material  in  terms  of 
a few  simple  tensile  tests  is  obvious.  However,  several  potential  limita- 
tions of  the  tensile  test,  are  widely  recognized.  These  include  (1)  limited 
strain  range  due  to  necking  instabilities,  (2)  limited  strain  rate  range 
below  most  forming  rates,  and  (3)  lack  of  information  on  stress  state  depend- 
ence. In  general,  these  limitations  have  not  had  a negative  impact  on  the 
testing  program,  and  this  is  discussed  in  detail  where  appropriate  in  the 
report.  In  regard  to  (3)  above,  it  has  been  possible  to  compare  the  uniaxial 
test  results  with  some  biaxial  tests  obtained  in  related  programs. 

The  constitutive  relations  determined  can  be  used  to  calculate  theore- 
tical forming  limit  curves  (FLC).  The  FLC  has  been  calculated  based  upon 
Hill's  [3]  condition  for  localized  necking.  The  limitations  of  the  calcu- 
lation will  be  discussed. 

1.4  Relation  to  Other  Air  Force  Programs 

This  program  has  been  conducted  in  close  cooperation  with  two  other 
Air  Force  programs  concerning  sheet  metal  formability  and  sheet  forming 
processes. 

(1)  Work  on  the  analytical  determination  of  forming  limit  curves  was 
conducted  by  Dr.  Louis  J.  Teutonico  as  a visiting  scientist  and 
consultant  to  the  Air  Force  Materials  Laboratory.  The  forming 
limit  curve  analysis  used  here  was  developed  under  this  effort. 

Dr.  Teutonico' s work  has  continued  as  part  of  contract  F33615- 
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78-R-5025  entitled  "Research  to  Develop  Process  Models  for 
Producing  a Dual  Property  Titanium  Alloy  Compressor  Disk" 
with  Battelle  Columbus  Laboratories  as  prime  contractor. 

(2)  A program  entitled  "Mathematical  Modeling  of  Sheet  Metal  Forma- 
bility  Indices  and  Sheet  Metal  Forming  Processes"  is  being  con- 
ducted under  contract  F3361 5-77-C-5059  by  Battelle  Columbus  Lab- 
oratories as  part  of  the  Air  Force  ICAM  Program.  ALCOA  Labora- 
tories, University  of  Kentucky,  University  of  California  (Berke- 
ley), and  McDonnell  Aircraft  Company  are  subcontractors.  The 
tensile  test  matrix  performed  here  has  been  coordinated  with  this 
program.  Experimental  forming  limit  diagrams  (FLD's)  on  2024-0 
Av  (by  ALCOA)  and  on  T l - 6Av -4 V (by  UK)  are  available  for  compar- 
ison to  theoretical  predictions.  The  tensile  test  results  and 
FLD's  referenced  in  this  report  are  described  in  detail  in  the 
first  annual  report  for  the  Battelle  contract  [4].  Bulge  test 
results  on  2024-0  At’  by  ALCOA  completed  after  the  annual  report 
date  have  been  transmitted  to  the  author  as  a private  communi- 
cation [5]. 
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2.  CONSTITUTIVE  RELATIONS 


2 . 1 Strain  Hardening  Relationships 

In  most  metal  forming  calculations,  the  constitutive  relations  employed 
are  empirical  equations  which  relate  flow  stress  to  strain  or  strain  rate  hut 
usually  not  both.  When  strain  and  strain  rate  are  considered  together,  it  is 
often  in  terms  of  an  empirical  strain  hardening  law  with  a superposed  power 

• m 

law  creep  strain  rate  factor  k . 

The  most  ubiquitous  empirical  strain  hardening  relationship  is  the 
power  law  suggested  by  Holloman  [6], 

o = k (2-1) 

The  parameter  n is  the  strain  hardening  exponent  which  is  the  cornerstone 
of  many  metal  forming  calculations.  However,  Eq.  (2-1)  is  usually  an  over- 
simplification in  that  the  strain  hardening  exponent  is  not  a constant  hut 
depends  upon  strain  and  strain  rate.  This  leads  to  a strain  hardening  law 
of  a different  form  than  Eq.  (2-1). 

Other  strain  hardening  laws  have  been  proposed.  Those  discussed  most 
often  are  by  Ludwik[7], 

a = oQ  + k c\  (2-2) 

by  Swift  [8], 

-■M.  VO".  (?-« 

and  by  Voce  [9,  10], 

o ::  n - (o  - i>0)  exp  (-i/n).  (2-4) 

The  main  advantage  of  the  Ludwik  and  Swift  laws  over  the  power  law  is 
simply  that  they  Include  three  parameters  rather  than  two  and,  hence, 
may  account  for  a wider  range  of  strain  hardening  behavior  such  as  a 
finite  yield  stress  or  effects  of  a prestrain  i . However,  It  should  be 
noted  that  the  Swift  equation  will  appear  in  the  analysis  in  this  report 
from  a somewhat  different  perspective. 

The  Voce  equation  Includes  not  only  a finite  yield  stress  ciQ , but 
also  a saturation  stress  o . for  this  reason  among  others,  Kncks  [11] 
employed  the  Voce  equation  to  describe  the  strain  rate  dependent  strain 


hardening  behavior  of  high  purity  aluminum  and  typo  304  stainless  stool. 

In  order  to  fit  the  Voce  equation,  he  neglected  the  high  strain  portion 
of  the  tensile  data.  The  strain  hardening  behavior  to  bo  reported  here, 
particularly  for  2024-0  At,  does  not  support  the  concept  of  a saturation 
stress  and  suggests  that  Kocks'  approximation  in  neglecting  the  high  strain 
tensile  data  may  not  have  been  justified. 

An  important  feature  of  our  analysis  of  strain  hardening  behavior 
will  be  an  emphasis  on  the  strain  hardening  derivative  and,  in  particular, 
on  log-log  plots  of  the  strain  hardening  coefficient  a vs.  o.  A similar 
approach  lias  been  suggested  by  Reed-Hill  et  al . [I?]  who  plot  log  (de/di  ) 
vs.  log  e.  For  nickel  270  stress-strain  curves  obtained  at  300  K , they 
observe  stages  where  log  (de/di ) decreases  linearly  with  log  e which  is 
similar  to  the  results  to  be  reported  here. 

Our  experimental  results  show  clearly  that  the  strain  hardening 
coefficient  > and  strain  rate  sensitivity  v depend,  in  general,  on  both 
strain  and  strain  rate.  Hence,  strain  rate  independent  strain  hardening 
relationships  are  not  sufficient  to  model  the  plastic  behavior.  A more 
general  approach  is  required  to  guide  the  development,  of  constitutive 
relations . 

2.2  Hart's  Phenomenological  Approach 

The  overall  requirements  for  material  constitutive  relations  appli- 
cable to  sheet  metal  formability  have  led  us  to  the  phenomenological 
approach  developed  by  L.W.  Hart  [13]  and  recently  summarized  by  Hart  ct  al. 
[14],  This  is  a flow  theory  in  which  the  flow  stress  is  expressed  as  a 
function  of  strain  rate,  absolute  temperature,  and  one  or  more  explicitly 
history-dependent  parameters  which  characterize  the  current  structure  of 
the  material  and  evolve  with  continuing  deformation,  These  structure 
parameters  partially  determine  the  current  mechanical  properties  and  can 
be  used,  for  example,  to  describe  lot-to-lot  variations.  Strain  hardening 
can  be  described  in  terms  of  the  evolution  of  the  structure  parameters 
with  accumulated  plastic  strain.  Another  important  feature  of  this  approach 
is  that  the  two  most  important  formability  indices,  the  strain  hardening 
coefficient  > and  the  strain  rate  sensitivity  v,  appear  as  state  functions 
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such  that  their  dependence  on  stress,  strain  rate,  and  temperature  can 
be  specified.  This  Implies  that  the  Incremental  stress-strain  relationship 

df.no  = id.  -*  v dtn*  ( ;’-‘j ) 

can  be  integrated  for  any  known  deformation  path. 

An  analytical  form  for  the  rate  dependence  of  the  flow  stress  has 
recently  been  proposed  by  Hart  [15],  namely 

O = o'*  exp  [ -(  7 ) ] 4 i’0  > 1/M.  (?-6) 

I 

This  equation  represents  non-elastic  grain  matrix  deformation  only  with 
the  first  term  describing  thermally-activated  flow  and  the  second  term 
describing  dislocation  glide.  The  parameters  \ and  M are  material  constants 

•k  , ★ 

and  .',i.  , and  o represent  the  current  mechanical  state.  The  ovoluntioiuiry 

* # * 

parameters  <>  and  . are  related  so  Fg.  (2-6)  includes  two  independent  param- 
eters. The  equation  should  be  applicable  at  all  temperatures  below  about 
0.4  times  the  melting  temperature  above  which  grain-boundary  sliding  must 
also  be  considered. 

The  problem  of  developing  an  analytical  form  to  describe  strain  hard- 
ening In  terms  of  the  evolution  of  the  history-dependent  parameters  in  fg. 

( 2-6 ) is  not  as  well  developed.  It  is  expected  that  the  extensive  experi- 
mental data  on  strain  hardening  vs.  temperature  and  strain  rate  obtained 
in  this  program  will  contribute  to  a solution  to  this  problem. 
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3.  FORMING  LIMITS 


3. 1 Experimental  Forming  Limit  Diagrams 

The  forming  limit  diagram  (FLO)  approach  to  assessing  limiting  strains 
to  failure  in  a sheet  metal  process  was  introduced  by  S.P.  Keeler  [16]  in 
1965.  The  FLO  presents  a locus  of  principal  strains  (e^ , e^)  in  the  plane 
of  a sheet  beyond  which  failure  is  expected.  By  definition,  e^  is  the 
largest  principal  strain,  and  we  are  concerned  with  failure  related  to 
thinning  (e^  < 0).  Hence,  the  FLD  covers  strain  states  from  pure  shear 
(true  strain  =-e^)  through  plane  strain  (e2  3 0)  to  balanced  biaxial 
tension  (e-j  = e^). 

The  FLD  is  determined  experimentally  by  laboratory  tests  on  gridded 
blanks  which  produce  the  appropriate  range  of  principal  strain  ratios  e^/e-j . 
Keeler's  original  work  [16]  employed  biaxial  stretching  experiments  to  pro- 
duce tension-tension  strain  states  to  determine  the  right-hand  side  (e2  > 0) 
of  the  diagram.  Goodwin  [17]  extended  this  work  to  the  left-hand  side 
(e^  < 0)  using  various  types  of  cup  and  tension  tests  for  tension-compres- 
sion strain  states.  The  FLD  is  defined  in  terms  of  the  major  and  minor 
strains  measured  on  a grid  of  deformed  circles  at  the  onset  of  visible, 
localized  necking.  It  is  drawn  to  fall  below  the  strains  in  necked  and 
fracture-affected  zones  and  above  the  strains  found  just  outside  these 
zones.  The  original  Keeler-Goodwi n FLD  was  presented  as  a narrow  band 
which  helped  account  for  the  experimental  uncertainty  in  identifying 
the  forming  limits. 

The  techniques  currently  being  used  to  generate  experimental  FLD's 
can  be  divided  into  two  groups,  in-plane  and  out-of-plane.  For  in-plane 
techniques,  the  sheet  is  not  bent,  and  there  are  usually  no  surface 
(frictional)  forces  cr  normal  pressures.  These  techniques  involve  modi- 
fied tensile  specimens  and  punch-stretching  of  a reduced  thickness  patch 
[18].  The  more  standard  out-of-plane  techniques  usually  follow  the  modi- 
fied punch-stretch  test  routine  suggested  by  Hecker  [19],  In  these  tests, 
a gridded  sheet  is  clamped  securely  at  the  periphery  and  stretched  to 
failure  over  a hemispherical  punch.  Failures  are  generated  over  a range 
of  strain  ratios  e^/e^  by  increasing  the  lubrication  on  the  positive  side 
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(e,,  > 0)  and  decreasing  the  blank  width  on  the  negative  side  (e0  <•  0). 

The  deformation  modes  in  the  out-of-plane  tests  of  course  correspond 
more  closely  to  actual  sheet  metal  processes. 

While  there  are  some  cases  for  which  in-plane  and  out-of-plane  FLD's 
agree,  it  is  generally  found  that  the  out-of-plane  limits  are  appreciably 
higher  than  the  in-plane  limits.  This  has  been  discussed  by  Ghosh  and 
Hecker  [20]  and,  more  recently,  by  Ghosh  [21].  The  differences  are  attri- 
buted to  variations  in  strain  path  and  in  the  process  of  strain  localization 
which  is  much  slower  for  out-of-plane  tests. 

3 . 2 Theoretical  Forming  Limit  Curves 

Although  the  FLO  is  a basic  and  useful  way  to  assess  formability,  it 
would  he  difficult  to  determine  the  explicit  dependence  of  limit  strains 
on  material  properties  or  even  identify  which  properties  had  the  most  direct 
affect.  Hence,  an  important  complimentary  technique  to  the  Fin  is  the 
theoretical  forming  limit  curve  (FLC).  The  FLC  represents  'cific  failure 
criterion  and,  hence,  a specific  set  of  relevant  property 

Plastic  Instability  Criteria.  The  most  important  failure  criteria  for 
sheet  forming  describe  some  mode  of  plastic  instability.  Swift  [8]  dis- 
cussed the  onset  of  diffuse  necking  and  showed  that  it  can  begin  when  the 
strain  hardening  coefficient  > reduces  to  a critical  value  which  depends  on 
the  strain  ratio  For  an  isotropic  sheet  deformed  by  uniaxial 

stress  i 2^  = - 1/2  and  the  critical  value  is  y * 1.  However,  the  onset 
of  diffuse  necking  is  not  readily  observable,  and,  hence,  the  Swift  cri- 
terion is  nut  a relevant  failure  criterion  for  sheet  metal. 

Plastic  instability  leads  to  failure  only  when  a localized  neck 
forms.  Hill  [3]  presented  the  .ondition  for  the  onset  of  a localized  neck 
in  a thin  sheet.  Again,  the  condition  is  stated  in  terms  of  a critical  ■> . 

For  an  isotropic  sheet  deformed  by  uniaxial  stress  it  is  ^ = 1/2.  However, 
the  Hill  condition  applies  only  to  the  left  hand  side  of  the  FI.C  for  which 
«-  0.  Also,  it  does  not  accourt  for  strain  rate  dependence. 

To  account  for  localization  of  deformation  for  strain  states  with 
i 2 > Of  two  approaches  have  been  taken.  The  first,  due  to  Marciniak  and 
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Kuczynski  [2?]  ( M- K ) and  extended  by  Hutchinson  and  Neele  [2:1],  postulates 
ilie  existence  of  an  initial  non-homogeneity  in  the  form  of  a thickness 
groove  across  the  sheet.  As  the  sheet  deforms,  the  strain  becomes  local- 
ized to  this  groove.  The  second  approach,  first  proposed  by  Storen  and 
Rice  [24]  and  later  modified  by  Hutchinson  and  Neale  [23],  incorporates 
a deformation  theory  of  plasticity  into  a classical  bifurcation  analysis. 

The  bifurcation  mode  corresponds  to  localized  deformation  in  a narrow  band, 
as  in  Hill's  analysis.  Only  the  first  approach  will  he  discussed  further. 

The  M-K  analysis  was  based  upon  a von  Mises  yield  function  incorpo- 
rating normal  anisotropy  and  Swift's  strain  hardening  relation,  Tq.  (2-3). 
Strain  rate  effects  were  not  included.  The  M-k  approach  was  also  investi- 
gated by  Sowerby  and  Duncan  [25]  who  discussed  the  dependence  of  forming 
limits  on  material  properties  and  by  Lee  ana  icobayasni  [2]  who  extended  the 
calculation  to  the  left-hand  side  (i0  •-  0)  and  discussed  planar  anisotropy 

L 

and  strain  path  dependence. 

Strain  rate  effects  were  first  introduced  by  Marciniak,  Kuczynski. 

and  Pokora  [26]  (M-K-P)  uid  discussed  more  recently  by  Hutchinson  and 

[23],  The  strain  rate  - pendence  is  introduced  through  the  constitut 
, m 

relation  by  including  a pu*or  law  creep  facto.  . . It  is  shown  that  quite 

small  values  of  the  strain  rate  sensitivity  (m  ; 0 01 ) can  haze  a large 
positive  effect  on  the  forming  limits.  The  strain  rate  dependent  calcu- 
lations have  been  developed  by  Ghosh  [21]  who  found  good  agreement  between 
calculated  and  experiim  _al  forming  limits  for  A-K  steel.  The  rate  sensi- 
tivity provided  nearly  a 50  percent  enhancement  in  the  calculated  plane- 
strain  forming  limit. 

There  are  two  problems  with  the  M-K  and  M-K-P  analyses  which  deserve 
mention.  One  is  the  requirement  for  an  initial  inhomogeneity.  This  ts 
usually  expressed  in  terms  of  a thickness  (f  = At/t)  or  strength  (f  ~ Ak/k, 
k from  Fq.  (2-3)) defect  or  combination  of  both.  The  calculated  forming 
limits  are  quite  sensitive  to  the  assumed  size  of  the  initial  inhomogeneity. 
Ghosh  [21] finds  agreement  between  theory  and  experiment  for  f r 0.002  - 0.010 
which  he  justifies  as  being  reasonable  values.  Still  f is  not  directly 
measureable  and  must  be  viewed  at  present  as  an  adjustable  parameter.  The 
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other  problem  involves  the  extension  of'  the  analysis  to  the  left-hand  side 
of  the  FIX.  On  the  right-hand  side,  the  groove  should  be  taken  perpendic- 
ular to  the  avis  of  the  largest  principal  strain.  For  the  left-hand  side 
this  is  not  the  case.  Hutchinson  and  Neale  [03]  discuss  the  problem  of 
optimizing  the  angle  between  the  groove  and  the  principle  strain  direction 
so  as  to  obtain  the  lowest  forming  limit.  For  t he  strain  rate  independent 
case,  it  is  probably  a good  approximation  to  take  the  groove  (defect)  at 
the  Hill  [37]  angle  though  for  the  strain  rate  dependent  case  this  may 
not  be  justified  [?$].  Since  these  problems  have  not  been  resolved,  the 
FLC  calculations  in  this  report,  will  he  based  on  the  Hill  condition.  Un- 
fortunately, this  means  that  strain  rate  effects  are  not  included. 


Hill's  Condition  for  Localized  Necking.  The  derivation  of  the  FLC 
based  on  Hill's  theory  of  localized  necking  is  described  in  reference 
[!].  It  will  be  reviewed  briefly  here.  Hill's  condition  is  [3], 

(3-1) 


(£♦-)/  4 . 

A , .3  <?  X' 


1 

where  f is  the  yield  function  as  defined  in  Appendix  A.  Other  assumptions 
are  plane  stress,  normal  isotropy,  and  proportional  straining.  Using  Lgs. 
(A4,  AS,  A20)  of  Appendix  A,  it  can  be  shown  that  the  FLC  is  given  by 


1 1 


i ,,  - > i 


(3-3) 


The  right  hand  side  in  Eg.  (3-3)  can  be  determined  from  the  constitutive 
relation  as  a function  of  strain  ratio  (:—,/<  ^ . 

To  plot  the  FLC,  we  use  a parametric  form  determined  from  [gs. 

(A4,  AS).  For  the  Swift  strain  hardening  relation,  Lg.  (3-3),  the 
effective  strain  . is  related  to  the  critical  > by 


X = n / (.o«  D. 


(3-3) 
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4.  EXPERIMENTAL  TEST  PROGRAM 


4 . 1 Material s and  Tost  Specimens 

The  2024-0  aluminum  and  titanium  -6Av-4v  sheet  materials  used  in  this 
project  were  supplied  through  Pattelle  Columbus  Laboratories  from  lots 
obtained  for  their  sheet  metal  forming  project  (see  sec.  1.4).  Both  sheet 
materials  wore  0.030  in  (1.27  mm)  thick. 

The  2024-0  A;  sheet  was  fabricated  by  Reynolds  Metals  Company  and 
supplied  according  to  ML  specif icat ion  Mil -QQA-230/4 . It  was  not  clad. 

The  grain  size,  determined  by  ALCOA  141,  was  a relatively  large  ASTM  3. -3. 

The  T i - 1'  A v - 4 V sheet  was  fabricated  by  RMI  Company  and  supplied  accord- 
ing to  MIL  specification  MIL-1-9046,  type  111,  composition  C in  the  annealed 
condition  (1450'T  for  13  min.  + A.C.).  Metal lographic  examination  in  our 
laboratory,  at  the  University  o.  Kentucky  [ 4],  and  at  Pattelle  [ 4] 
indicated  that: 

1.  The  microstructure  is  characterized  by  a fine  grain  size  of 
2-3  cm. 

2.  The  material  was  not  fully  (reoystai  1 ized)  annealed.  This  is 
consistent  with  the  MIL  specification  and  industry  practice. 

3.  Some  plastic  work,  as  evidenced  by  apparent  grain  elongation  along 
the  rolling  direction,  remained  in  the  material . 

The  crystal lographic  texture  of  the  Ti-oA<-4V  sheet  was  also  determined 
using  X-ray  polo  figures  by  Boeing  Technology  Services,  Seattle,  Washington, 
and  reported  bv  Pattelle  [ 4 1 . Interpretation  of  the  (0002)  pole  figure 
indicated  that  three  textural  components  were  dominant.  (11  basal  poles 
aligned  with  the  transverse  direction,  (2)  basal  poles  rotated  approximately 
131'  from  the  sheet  normal  toward  the  rolling  direction,  and  (3)  basal  poles 
aligned  with  the  rolling  direction.  Neither  component  was  particularly 
strong  with  intensitie  ot  appi\  ..liatelv  .'x  random..  Tims  is  a coi.Jiion 
conmerical  li-bA.-4V  tenure  aiul  leads  to  r-values  1. 

The  sheet  materials  were  cut  into  test  specimen  blanks  which  were 
machined  into  pin-loaded  tensile  specimens  ot  1 . t->0  in  ^3G.l  mm)  guago 
length.  The  specimen  ursign  is  shewn  in  lig.  1.  All  dimensional  ratios 
correspond  to  the  pm-leaded  sheet  specimen  in  AS1M  sooc i * ica t ion  I-;''. 
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NOTES: 

1.  Width  of  reduced  section  must  be  accurate  to  ±0.002  and 
uniform  to  +0.0005  within  1.50  in  reduced  length. 

2.  Holes  must  be  on  of  reduced  section  to  ±0.0001. 


figure  i.  pin-loaded  sheet  tensile  specimen 
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Before  cutting  the  sheets,  maps  and  specimen  codes  were  prepared  so  that  the 
location  of  each  specimen  in  the  sheet  could  he  subsequently  determined. 
Tensile  specimens  were  cut  from  the  sheet  in  three  different  orientations, 
at  0°  (L),  90''  (T),  and  45''  (F)  to  the  rollinq  direction. 

It  is  clear  from  the  results  in  Sec.  6 of  this  report  that  the 
T i - 6 A v - 4 V sheet  material  has  very  limited  strain  hardening  capacity.  It 
is  thought  that  this  is  due  to  the  fact  that  it  was  only  partially  annealed 
and  retains  a structure  influenced  by  the  retention  of  previous  plastic 
work.  Hence,  for  the  study  of  the  strain  hardening  behavior  of  Ti-6Av-4V, 
some  experiments  have  been  performed  on  button-head  tensile  specimens 
machined  from  fully  annealed  T i -6A c - 4 V round  bar  stock  available  from  a 
previous  study  [L’91.  This  material  had  been  supplied  according  to  MIL 
specification  Ml L - T -9047  and  had  been  annealed  at  14B0T  for  ?.  hours  plus 
air  cooled. 

A,  2 Tens 1 1 e Test  Procedures 

Both  constant  extension  rate  and  load  relaxation  tests  were  carried  out 
using  a ?['  kN  Instron  HITT  screw-driven  tester.  The  crosshead  speed  could 
be  varied,  essentially  continuously,  from  BOO  mn/min  to  below  BxlO"'  nin/min. 
This  allowed  initial  strain  rates  of  10* 1 to  10"'  sec"'  to  be  obtained. 

For  load  relaxation  tests,  the  loading  strain  rate  was  chosen  to  be 
10*-’  sec-1.  Load  relaxation  data  was  recorded  after  stopping  the  instron 
crosshead  at  a predetermined  extension  as  continued  plastic  strain  in  the 
specimen  replaced  elastic  strain  in  both  the  specimen  and  load  train. 

For  tests  above  room  temperature,  specimen  temperatures  were  main- 
tained using  a 3-2one  split  tube  furnace  with  a special  high  stability 
Furotherm  three  mode  controller.  For  load  relaxation  tests  and  for  tensile 
tests  at  a strain  rate  of  10*’'  sec*1  which  lasted  up  to  10  hours,  tempo  a - 
tures  were  controlled  to  better  than  *0.i>"C. 

The  determination  of  plastic  strain  ratios  was  carried  out  by  two 
different  methods.  (1)  Direct  width  (w)  and  thickness  (t)  measurements 
were  made  after  interrupted  tests  and  compared  to  initial  values.  In 
this  case. 


fPnnf|l||f1IF'» 


r = 


Lw 


(2)  For  room  temperature  tests,  an  Instron  transverse  strain  sensor  was  used 
to  record  the  decrease  in  thickness  continuously.  In  this  case,  assuming 
constant  volume  during  plastic  deformation, 


The  factor  de  /dc,,  can  be  determined  from  the  slope  of  the  transverse  strain 
t t 

sensor  output. 


4.3  Automatic 


For  both  constant  extension  rate  and  load  relaxation  tests,  load-time 
data  was  obtained  digitally  using  a Hewlett-Packard  (HP  3052A)  automatic 
data  acquisition  system.  Digital  data  acquisition  is  a necessary  component 
of  the  testing  program  in  order  to  provide  sufficiently  rapid  data  acqui- 
sition at  the  highest  strain  rates  and  sufficient  load  resolution  at  the 
lowest  strain  rates  in  load  relaxation  tests  and  near  the  maximum  load  in 
constant  rate  tests.  The  resolution  is  necessary  to  be  able  to  calculate 
accurate  derivatives  needed  to  investigate  strain  hardening  and  load 
relaxation . 

The  HP  3052A  data  acquisition  system  consists  of  a scanner,  digital 
voltmeter  (DVM),  timing  generator,  and  plotter  controlled  by  a calculator 
through  a coninon  interface  bus.  A block  diagram  of  the  system  is  shown  in 
Fig.  2.  Data  and  programs  can  be  displayed  on  a line  printer  or  stored  on 
a magnetic  tape  cassette  both  of  which  are  components  of  the  calculator/ 
controller.  Output  can  also  bp  sent  to  a teletype  (TTY)  through  a serial 
interface  (RON)  module.  The  system,  as  assembled,  has  15,036  bytes  of  read/ 
write  memory,  and  problems  due  to  limited  storage  can  usually  be  overcome  by 
using  the  tape  cassettes. 

A program  (DVM  RUN)  has  been  developed  which  enables  the  calculator  to 
control  the  scanner  and  DVM  reading  sequence.  The  channels  scanned,  reading 
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BLOCK  DIAGRAM  Of  THE  AUTOMATIC  DATA  ACQUISITION  SYSTEt 


intervals,  and  output  devices,  including  changes  in  these  instrument  settings, 
can  be  preselected  for  a particular  experiment.  A flow  chart  for  the  DVM  RUN 
program  is  shown  in  Fig.  3.  The  main  program  accepts  a data  sequence  matrix, 
triggers  the  timing  generator,  sets  up  data  storage  and  output.  Each  timing 
generator  pulse  initiates  an  interrupt  routine  (CSI)  which  either  stores  a 
DVM  reading  or  changes  the  instrument  settings.  The  maximum  reading  rate 
depends  somewhat  on  the  data  sequence  but  is  approximately  15  sec . Since 
the  system  has  a buffered  output  capability  (limited  by  the  15k  core  size), 
the  reading  rate  is  not  affected  by  the  speed  of  the  output  device. 

The  DVM  readings  are  stored  in  a data  string  (D$).  Interpretation  of  the 
DVM  readings  by  channel  (load,  strain,  etc.)  and  time  (after  t=0)  is  accom- 
plished by  a second  program,  DVM  DATA,  which  identifies  the  DVM  readings  (Q$) 
according  to  the  data  sequence  matrix  (S$)  and  stores  voltage-time  data 
pairs  by  channel  in  (C$).  The  flow  chart  for  DVM  DATA  is  shown  in  Fig.  4. 

At  present,  the  program  DVM  DATA  assumes  that  up  to  two  channels  will  be 
used,  but.  this  could  easily  be  increased.  In  addition,  some  of  the  informa- 
tion in  (S$)  is  also  stored  in  (D$),  redundancy  that  was  convenient  during 
program  development. 

Listings  cf  the  programs  DVM  RUN  and  DVM  DATA  are  pr  sented  in  Appendix 
B.  The  language  is  a modified  form  of  Basic. 
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5.  RESULTS  FOR  2024-0  ALUMINUM 


5 . 1 Constant  Extension  Rate  Tests 

Stress-strain  curves  and  tensile  test  parameters.  Constant  extension 
rate  tensile  tests  were  performed  on  the  2024-0  aluminum  test  specimens  at 
temperatures  of  25°,  100°,  163°  and  225°C.  The  163°C  (325°F)  temperature 
was  selected  to  represent  a typical  warm  forming  temperature  for  this  alloy. 

At  25°C,  tensile  tests  were  performed  on  0°(L),  90°(T),  and  45°(F) 
specimens  at  strain  rates  of  10"5  to  10"1  sec'1.  For  these  tests,  the 
Instron  transverse  strain  sensor  was  used  to  monitor  the  decrease  in 
specimen  thickness.  A typical  Instron  chat t record,  obtained  at  a strain 
rate  of  10‘3  sec'1,  is  shown  in  Fig.  5.  At  strain  rates  of  10‘:'  sec"1  and 
below,  the  25°C  load-elongation  curves  showed  strong  load  serrations 
(Portevin-LeChatel ier  effect),  and  the  thickness  - elongation  curves  showed 
highly  non-uniform  flow  characteristic  of  Luders  band  propagation.  The 
load  serrations  only  began  after  a critical  strain  which  increased  with  in- 
creasing strain  rate.  At  10*1  sec"1,  the  critical  strain  was  high  enough 
that  only  one  or  two  serrations  were  ohserved.  While  these  effects  are 
well-known  in  aluminum  copper  alloys  [30]  and  have  recently  been  reported 
for  2024  aluminum  in  particular  [31],  use  of  the  thickness  strain  trans- 
ducer appears  to  have  produced  a unique  and  exceptionally  clear  picture  of 
the  band  propagation  phenomenon.  In  spite  of  the  localized  deformation,  the 
material  appears  to  strain-harden  in  a normal  manner,  and,  for  the  purpose 
of  studying  formability,  we  can  consider  an  averaged,  or  smoothed,  load- 
eiongation  curve. 

The  load -elongation  curves  have  been  used  to  compute  a number  of  stan- 
dard tensile  test  parameters.  These  are  summarized  for  the  2'a°C  tests  in 
Table  1.  In  this  and  subsequent  tables,  <<  is  the  engineering  stress  at 
0.2  percent  non-elastic  strain,  is  the  engineering  stress  at  maximum 

load,  i is  the  true  non-elastic  strain  at  maximum  load,  and  i ^ is  the 
engineering  strain  at  fracture  (1,50  in  guaye  length). 
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FIGURE  5.  LOAD  vs.  ELONGATION  AND  CHANGE  IN  THICKNESS  vs 
ELONGATION  FOR  2024-0  ALUMINUM  AT  25°C  AND  10“ 


TABLE  1.  TENSILE  TEST  PARAMETERS  FOR 


2024-0  ALUMINUM  AT  25"C 


Specimen 

Rate  (sec-1) 

0 (MPa) 

°UTS 

‘ur,  if 

tot^ 

19L 

lO'1 

90.5 

181  .5 

0 . 1 38 

17.7 

17L 

10‘p 

91.5 

182.1 

0.153 

20.4 

15L 

10- 3 

88.6 

182.7 

0.147 

18.3 

16L 

io-“ 

92.0 

186.9 

0.151 

19.2 

18L 

10's 

91.1 

191  .0 

0.160 

20.5 

20T 

io-> 

91.7 

171  .0 

0.141 

17.4 

13T 

10" 

89.7 

172.4 

0.157 

18.8 

14T 

io- '> 

92.3 

179.1 

0.149 

18.7 

12F 

10'1 

89.9 

171.3 

0 . 1 50 

18.9 

1 OF 

10*1 

85.4 

173.1 

0.157 

19.4 

1 IF 

10-,J 

90.3 

180.9 

0.156 

20.9 

In  Table  1,  the  stress  appears  to  be  Independent  of  strain  rate 
for  each  orientation.  The  small  variations  obtained  may  be  due  to  problems 
in  correcting  for  the  effect  of  load-train  compliance.  The  i ^ 1 

values  show  some  scatter  but  also  seem  independent  of  strain  rate.  How- 
ever, the  values  appear  to  show  a snail,  yet  monotonic  and  probably 
significant  decrease  with  increasing  strain  rate.  This  Is  believed  to  be 
simply  another  feature  of  the  mechanism  which  produced  the  localized 
plastic  deformation  illustrated  in  Fig.  5.  In  effect,  the  flow  stress  is 
determined  by  the  strain  rate  characteristic  of  the  band  propagation,  which 
is  different  from  and  not  simply  related  to,  the  rate  imposed  by  the 
testing  machine. 

For  tests  at  any  one  strain  rate,  the  values  for  the  1 and  f 
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specimens  are  nearly  equal  and  significantly  less  than  those  for  the  L 
orientation.  This  will  be  discussed  In  more  detail  below  in  terms  of 
plastic  anisotropy.  The  0^^.^  values  reported  here  are  in  excellent 
agreement  with  those  obtained  by  ALCOA  [4]  (L : 182  MPa,  T : 1 72  MPa)  at  a 
strain  rate  of  5x10‘3  sec"1  on  specimens  obtained  from  the  same  Reynolds 
2024-0  aluminum  lot. 

At  1 63°C , tensile  tests  were  also  performed  on  L,  T,  and  F specimens 
at  strain  rates  of  10-5  to  10'1  sec"1.  The  load-elongation  curves  were 
smooth  with  no  visible  Indication  of  localized  flow.  The  variation  of  ten- 
sile test  parameters  with  strain  rate  and  orientation  is  summarized  in 
Table  2. 


TABLE  2.  TENSILE  TEST  PARAMETERS  FOR 
2024-0  ALUMINUM  AT  163°C 


Specimen 

Rate  (sec*1) 

oy  (MPa) 

°UTS  {MPa) 

k unit 

‘ tot  ^ ^ 

51 L 

10-1 

87.2 

164.0 

0.120 

16.2 

4 CL 

10”2 

83.7 

162.6 

0.127 

21.6 

45L 

10* 3 

82.6 

155.0 

0.127 

33.3 

49L 

10-** 

81  .2 

138.9 

0.112 

30.6 

501 

10'b 

81.4 

131  .4 

0.121 

21  .9 

24T 

10'1 

79.7 

157.3 

0.122 

16.2 

23T 

10- 2 

81.4 

156.3 

0. 1 35 

21.0 

19T 

10*3 

83.9 

152.3 

0.130 

32.8 

18T 

10’u 

80.6 

136.3 

0.134 

34.8 

40F 

10-' 

81.5 

154.9 

0.124 

17.2 

13F 

10- 2 

80. B 

154.5 

0.146 

24.7 

25F 

10- 3 

81.4 

147.0 

0.137 

34.4 

28F 

10-'* 

79.0 

133.9 

0.135 

35.0 

Again  and  appear  to  be  effectively  independent  of  strain  rate. 

However,  the  values  now  indicate  a clearly  positive  strain  rate 
sensitivity,  and  the  ^ values  are  strongly  strain  rate  dependent  with  a 
maximum  near  10'3  or  10'**  sec-1.  It  has  been  shown  [32]  that  the  post- 
uniform elongation  (t^ot  - c in  a tensile  test  is  determined  pre- 

dominantly by  the  strain  rate  sensitivity.  Hence,  the  variation  in  r 
in  Table  2 indicates  a rate  dependent  strain  rate  sensitivity  that  is 
largest  near  10" 3 or  lO"4  sec"1. 

The  variation  in  tensile  properties  with  strain  rate  at  163°C  is  also 
illustrated  in  Fig.  6 by  the  true  stress  - true  strain  curves  (up  to  e ,.^) 
for  the  five  L specimens  of  Table  2.  This  figure  also  suggests  that  the 
influence  of  strain  rate  is  largest  between  10~3  and  10" 14  sec"1. 

lensile  tests  at  100°C  and  225°C  were  conducted  only  on  L specimens  at 
a strain  rate  of  10";  sec"1.  The  true  stress-true  strain  curves  (up  to 
t if)  for  the  L specimens  at  10"-’  sec"1  for  all  four  test  temperatures  are 
shown  in  Mg.  /,  and  the  tensiie  properties  are  summarized  in  Table  3. 

The  o , and  values  all  decrease  with  increasing  temperature. 

The  large  i ^ value  at  225°C  is  an  indication  of  a large  strain  rate 
sensitivity  for  that  test. 


TABLE  3.  TENSILE  TEST  PARAMETERS  FOR 
2024-0  ALUMINUM  A1  10":  sec'1 


Specimen  Temperature  ("C)  o (MPa) 


‘ un  i f 


1 62 . 6 
1 03 . 1 


0.153 

0.151 

0.127 

0.080 


202H-O  RL  "L* 
IS3  C 


00 

FIGURE  6 


0.00  0(2 
TRUE  PLASTIC  STRAIN 

TENSILE  TRUE  STRESS  - TRUE  STRAIN  CURVES  FOR  2024-0 

aluminum  at  i63°c  vs.  s train  rate 


Plastir  Anisotropy.  The  plastic  anisotropy  of  the  aluminum  sheet  Is 
characterized  by  the  plastic  strain  rations  r = c^/c^  measured  for  specimens 
cut  at  0°(L),  90°(r),  and  45°(F)  to  the  rolling  direction.  The  r-values 
were  determined  from  direct  width  and  thickness  measurements  according  to 
Eq.  (4-1).  Measurements  were  made  during  interrupted  tests  conducted  at 
25°C  and  a strain  rate  of  10' 3 sec'1  and  at  strains  ot  approximately  2,  3,  5, 
and  10  percent.  No  significant  variation  with  strain  was  observed.  The 
measured  values  of  rQ,  r and  are  given  in  Table  4 along  with  values 
determined  by  Battelle  (at  lO'2  sec-1)  [4]  and  ALCOA  (at  5xl0'3  sec'1)  [4], 

TABLE  4.  R-VALUES  FOR  2024-0  ALUMINUM 
AT  25°C 


Orientation 

This 

Investigation 

— it 

(at  10  J sec  ) 

Battelle  [ 4 ] 

(at  10  ^ sec  ^ ) 

ALCOA  [ 4 1 
(at  5x10  ^ sec  ^ 

rO 

0.74 

0.45-0.62 

0.74 

r90 

0.64 

0.51-0.57 

0.56 

1 45 

0.77 

0.62-  .(.80 

“ “ 

_ 

It  should  also  have  been  possible  to  determine  the  r-values  from  the 
slopes  of  the  transverse  strain  sensor  output  (see  Fig.  5),  at  least  for 
strains  below  the  onset  of  unstable  flow.  However,  the  transverse  strain 
sensor  output  generally  indicated  larger  changes  in  thickness  and,  hence, 
smaller  r-values  than  those  measured  directly.  This  discrepancy  is  not 
i .Jerst.ood,  and  the  directly  measured  r-values  are  preferred.  However,  the 
linear  form  of  the  transverse  strain  sensor  output,  exclusive  of  the 
"staircase"  behavior  in  the  localized  flow  regime,  also  indicates  that  the 
r-values  do  not  depend  significantly  on  strain.  In  addition,  the  transverse 
strain  sensor  slopes  did  not  vary  noticeably  with  strain  rate  so  we  expect 
that  the  strain  rate  dependence  of  the  r-values  is  small. 

The  values  in  Table  1 show  that  the  tensile  flow  stress  is 

dependent  upon  the  orientation  of  the  test  specimen.  At  a strain  rate  of 
10~3  sec'1,  the  values  for  T and  F orientations  are  approximately  equal 
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but  are  about  6 percent  less  than  that  for  the  L orientation.  This  depend- 
ence can  be  discussed  In  terms  of  Hill's  [27]  theory  of  planar  plastic 
anisotropy.  The  relevant  yield  function  and  definition  of  effective  stress 
are  presented  in  Appendix  A.  In  particular,  expressions  for  the  ratios  of 
the  effective  stress  o to  the  tensile  stress  ay  in  terms  of  the  measured 
r-values  are  given  by  Eqs.  (A13-A15).  The  corresponding  relationships 
between  the  effective  strain  »'  and  tensile  strain  are  given  by  Eq.  (A17). 
The  same  formulation  can  be  used  to  relate  effective  stress  and  strain  to 
the  stress  and  strain  measured  in  a hydraulic  bulge  test.  These  ratios 

are  given  in  Lgs.  ( A 1 6 ) and  (A1S)  and  are  based  on  the  assumption  of  balanced 
biaxial  stress.  These  stress  ratios  and  strain  ratios  have  been  calculated 
for  the  measured  r-values  and  are  summarized  in  Table  5 . 


TABLE  5.  EFFECT  I VI  STRESS  - EFFECT  I V[  STRAIN 
CONVERSION  FACTORS  AT  25°C 


Test 

Or  ion tat  ion 

r / r 

Ten  si  le 

0 "(E) 

0.949 

1 .053 

Tons  i le 

90" (T) 

0.991 

1.009 

Tens  i le 

4 5 " ( F ) 

0.994 

1 . 006 

Bulge 

-- 

l .057 

0.946 

In  order  to  compare  the  stress-strain  curves  vs.  orientation,  the  effective 
stresses  at  effective  strains  of  0.05,  0.10,  and  0.15  are  given  in  Table  6 
for  the  three  10'-  sec'1  runs  of  Table  1.  In  addition,  bulge  test  results 
reported  by  ALCOA  [5]  on  specimens  obtained  from  the  same  Reynolds  2024-0 
lot  are  also  included.  The  bulge  test  stresses  are  average  values  from  four 
different  tests.  It  is  clear  from  Table  6 that  the  effective  stress-effective 
strain  curves  determined  for  the  three  tensile  orientations  are  effectively 
coincident . The  bulge  tost  results  are  not  as  close  and  differ  from  the  tensile 
results  by  about  10  percent.  Considering  the  approximate  manner  by  which  hydrau- 
lic pressure  and  change  in  thickness  are  converted  to  true  stress  - true  strain 
for  the  bulge  test  [5].  the  agreement  is  considered  sat isfactory.  However, 

• does  not  lend  support  to  the  anisotropy  correction  factors  as  the  tensile 


TABLE  6.  EFFECTIVE  STRESS  VS.  EFFECTIVE 
STRAIN  FOR  2024-0  ALUMINUM 
TENSILE  AND  BULGE  TESTS  AT  25°C. 


Test 

Specimen 

Effective  Stress  o at 

E = 

0.05 

0.10 

0.15 

Tensile 

15L  (0°) 

163.5 

186.4 

199.5 

Tensile 

13T  (90°) 

164.0 

185.2 

198.0 

Tensile 

1 OF  (45°) 

163.1 

185.1 

198.9 

Bulge 

ALCOA  [5] 

187.0 

206.2 

217.8 

results  do.  Similar  comparison  between  tensile  and  bulge  test  results  for 
aluminum  alloys  have  previously  been  reported  by  several  investigators  [33,  34]. 

It  is  concluded  that  Hill's  theory  of  planar  anisotropy  together  with 
r-values  determined  from  measured  strain  ratios  can  be  used  to  account  for  the 
plastic  anisotropy  of  2024-0  aluminum. 

Strain  Hardening  Behavior.  Since  the  true  stress  - true  strain  curves 
for  the  L,  T,  and  F orientations  at  25°C  can  be  brought  into  coincidence 
using  Hill's  theory  of  planar  anisotropy,  the  discussion  of  the  strain  harden- 
ing behavior  will  be  restricted  to  a single  (L)  orientation.  Referring  to 
Tables  1 and  2 it  is  seen  that  the  dependence  of  aUTS  on  orientation  is 
roughly  the  same  at  163°C  as  for  25°C,  with  aUTS  for  the  L orientation  greater 
than  that  for  T and  F.  So  it  will  be  assumed  that  the  r-values  measured  at 
25°C  can  also  be  used  at  163°C  and  the  tensile  data  at  both  temperatures  will 
be  reduced  to  effective  stress  - effective  strain  using  the  factors  in  Table  5. 

The  strain  hardening  behavior  will  be  investigated  in  terms  of  the  depend- 
ence of  the  derivative 

y = (o  en  o/D  e }•  (5-1) 

on  stress  at  various  strain  rates  c.  The  choice  of  variables,  y - y (o,e), 
is  based  on  Hart's  concept  [15]  that  these  are  state  variables  which  partially 
describe  the  current  structure  (and  future  response)  of  the  material.  Although 
the  experiments  are  at  constant  extension  rate  and  not  constant  strain  rate  as 
implied  in  the  definition  of  y , the  derivatives  of  the  stress-strain  curve  can 
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be  identified  as  y with  negligible  error.  For  a constant  extension 
rate  test  [35], 

d uno/de  = y - (3  Jtno /3  £n  e)  a y.-v.  (5-2) 

For  the  results  presented  here,  y>  1 and  v< 0 .1.  Hence,  the  strain  rate  sensi- 
tivity v makes  a negligible  contribution  to  the  total  derivative. 

Figure  8 shows  a plot  of  y vs  o for  the  2024-0  aluminum  test  run  at 
163°C  and  10"3  sec-1.  This  is  a smooth  curve  over  the  entire  range  of  data 
(equivalent  to  5xl0"3  < e < e .j^)  and  provides  no  obvious  clue  as  to  how 
to  model  the  strain  hardening.  It  is  more  instructive  to  examine  a plot  of 
log  y vs  log  a.  This  plot  is  shown  in  Fig.  9 for  the  same  test  as  Fig.  8. 
Regions  of  different  curvature  in  Fig.  8 now  appear  as  regions  of  different 
slope  in  Fig.  9.  The  most  distinctive  feature  is  the  linear  decrease  in 
log  y with  log  a at  large  stresses.  Reference  to  the  stress  strain  data  shows 
that  the  linear  portion  begins  at  a strain  of  approximately  e = 0.05.  Since 
we  are  interested  in  the  large  strain  behavior  for  predicting  formab'ility,  it 
is  sufficient  to  model  the  strain  hardening  for  e >0.05  phenomenologically 
by  assuming 

«,n  y = «.nA-B$.n  a 

where  A and  B are  constants  which  can  be  determined  by  linear  regression  from 
the  strain  hardening  data.  This  straight  line  is  shown  in  Fig.  9. 

Integrating  Eq.  (5-3)  we  obtain  an  expression  for  the  stress-strain  curve. 


5 = (AB)1/B  oB/AB  + e 1/8 


(5-4) 


where  o0  is  the  stress  at  e = 0.  Although  A and  B can  be  determined  from  the 
strain  hardening  data,  cr0  must  be  determined  using  the  stress-strain  data. 
While  there  are  several  ways  to  do  this,  we  have  used  a linear  regression 

— g 

analysis  of  a vs  e. 

It  is  to  be  noted  that  Eq.  (5-4)  is  the  same  form  as  that  originally  pro- 
posed by  Swift  [8], 

0 = k (eq  + e)n.  (5-5) 

However , In  our  application,  n = 1/B  is  determined  from  the  dependence  of  y 
on  stress,  and  c0  = a0B  / AB  need  not  be  interpreted  as  a prestrain.  In  fact, 
it  will  be  seen  that  for  2024-0  aluminum,  e0  will  come  out  negative.  This 
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simply  represents  the  faot  that  the  straight  line  fit  in  Fig.  8 extra- 
polates to  higher  derivatives  y at  low  stress  than  those  actually  measured. 
Hence,  the  fitted  Swift  equation  will  curve  below  the  experimental  points  at 
small  strains.  The  small  strain  regime  can  be  ignored  or  modeled  separately, 
as  required  by  the  application. 

Using  the  log  y vs.  log  a and  a vs.  e data,  the  Swift  equation  para- 
meters have  been  determined  vs.  strain  rate  for  the  five  L specimens  tested 
at  25°C  (see  Table  1)  and  at  163°C  (see  Table  2).  These  are  summarized 
in  Tables  7 and  8 respectively  along  with  a comparison  between  the  calcu- 
lated maximum  uniform  strain 

c .,  = n - e0  (5-6) 

umf  0 

and  the  measured  (effective)  value  eun^. 

The  success  of  this  fitting  procedure  is  illustrated  by  the  effective 
stress-strain  curve  and  Swift  equation  fit  for  the  163°C,  10  3 sec  1 test 
in  Fig.  10.  It  is  especially  significant  that  the  fitting  procedure  pre- 
sented here  produces  an  excellent  match  between  calculated  and  experimental 
slopes.  It  is  the  slope  of  the  stress  strain  curve  which  directly  deter- 
mines the  formability  when  failure  occurs  by  plastic  instability. 


TABLE  7.  SWIFT  EQUATION  PARAMETERS  FOR 
2024-0  ALUMINUM  AT  25°C. 


Specimen 

Rate  i 
(sec'  ) 

n 

i'0 

k(MPa) 

n-»o 

cunif 

19L 

10"1 

0.147 

-0.015 

266.8 

0.162 

0.145 

17L 

10'2 

0.134 

-0.023 

262.5 

0.157 

0.161 

15L 

O 

1 

co 

0.144 

-0.017 

266.7 

0.161 

0.155 

16L 

10"4 

0.116 

-0.032 

259.6 

0.148 

0.159 

18L 

-5 

10  3 

0.129 

-0.030 

273.1 

0.159 

0.168 

average 

values 

0.134 

-0.023 

265 . 7 

0.157 

0.158 

In  Table  7,  none  of  the  parameters  n,  to*  or  k for  the  25"C  tests  show 
a smooth  or  significant  variation  with  strain  rate.  The  same  is  true  for  the 
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uniform  strains.  Hence  it  is  concluded  that  the  high  strain,  strain  harden- 
ing behavior  of  2024-0  aluminum  at  25°C  is  effectively  independent  of  strain 
rate.  For  this  reason,  average  values  of  the  Swift  equation  parameters  and 
uniform  strains  have  been  included.  At  25°C,  the  constitutive  equation  based 
on  these  average  values  is 

o (MPa)  = 265-7  (-0.023  + T)0:134  (5.7) 

For  the  individual  runs,  and  particularly  for  the  average  values,  we  obtain 
excellent  agreement  between  calculated  and  observed  uniform  strains. 

TABIC  8.  SWIFT  EQUATION  PARAMFTERS  FOR 
2024-0  ALUMINUM  AT  1 63°C . 


Specimen 

Rate  , 
(sec'  ) 

n 

t ~ 

k(MI’a) 

n-,  0 

tuni  f 

51 L 

10' 1 

0.107 

-0.022 

223.6 

0.129 

0.126 

46L 

O 

i 

2-0 

0.109 

-0.021 

222 . 3 

0.130 

0.134 

45L 

10'3 

0.104 

-0.013 

onn  l 

0.122 

0.133 

49L 

o 

1 

0.098 

-0.016 

184.5 

0.114 

0.118 

50L 

10'5 

0.107 

-0.018 

178.3 

0.126 

0.127 

average  values 

0.105 

-0.019 

-- 

0.124 

0.128 

For  the  163°C  tests  in  Table  8,  neither  n,  i0,  or  the  uniform  strains 
show  a significant  variation  with  strain  rate  and  average  values  have  been 
presented.  However,  k increases  significantly  and  non-uni formly  with  increas- 
ing strain  rate.  Hence,  a simple  constitutive  equation  like  Lq.  (5-7)  will 
not  be  presented  for  the  163°C  strain  hardening  behavior.  The  strain  rate 
dependence  of  the  strain  hardening  behavior  is  further  illustrated  in  Fig.  11 
which  shows  (from  left  to  right)  the  log  > plots  for  strain  rates  of  10 
10~M,  10*',  and  10”‘'sec  _1.  The  10_l  sec*1  plot  has  been  omitted  since  it  is 
essentially  coincident  with  the  10*-"  sec'1  plot.  Figure  11  shows  that, 
although  the  n-values  are  constant  and  k-values  strain  rate  dependent  for 
the  high  strain  behavior,  essentially  the  opposite  is  the  case  for  the  low 
strain  behavior.  The  slight  upward  curvature  at  very  low  strains,  also 
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evident  in  Fig.  9,  represents  a transient  strain  which  appears  to  saturate 
by  e : 5x10  3.  The  strain  rate  dependence  of  the  strain  hardening  behavior  is 
largest  near  a strain  rate  of  10'3  or  10'1*  sec"1  which  is,  of  course,  con- 
sistent with  the  high  strain  tensile  test  results  presented  in  sec.  5.1.  In 
summary,  the  stress  and  strain  rate  dependence  of  > shown  tn  Fig.  11  indicates 
that  the  strain  hardening  behavior  of  2024-0  aluminum  at  163°C  cannot  be 
modeled  by  a strain  rate  independent  constitutive  equation.  The  explicit 
dependence  on  strain  rate  will  have  to  be  taken  into  account. 

Comparison  Between  Tensile  and  Bulge  Tests.  In  Table  6,  the  tensile 
effective  stress-effective  strain  data  is  compared  to  bulge  tests  results  on 
the  same  2024-0  lot  reported  by  ALCOA  [5].  This  comparison  is  made  within 
the  range  of  tensile  strains,  t < ’ ...  Since  one  important  objective  of 

the  development  of  constitutive  equations  is  the  ability  to  extrapolate  out- 
side the  range  of  tensile  data,  comparison  will  also  be  made  to  the  bulge 
tests  results  at  the  larger  strains  achievable  in  that  test. 

Figure  12  shows  a comparison  between  the  Swift  equation  fit  to  the  high 
strain  tensile  data  for  the  25°C,  10"  1 sec'1  test  and  several  bulge  test  data 
points.  The  comparison  is  similar  to  that  indicated  in  Table  6.  The  values 
of  effective  stress  determined  from  the  bulge  test  are  about  10  percent  above 
the  extrapolated  tensile  data.  The  strain  hardening  behavior  (slopes)  appear 
very  similar. 

The  strain  hardening  behaviors  are  compared  further  in  Fig.  13  which 
shows  log  Y vs.  log  a for  (a)  tensile  data  obtained  at  25°C,  10  1 sec  , 

(b)  the  average  Swift  equation  for  the  25°C  results  from  Eq.  (5-7),  and  (c) 
bulge  test  data.  It  is  seen  that  the  bulge  test  data  agree  well  with  the 
linear  extrapolation  of  the  decrease  in  log  > vs.  log  o based  on  the  high 
strain  tensi le  data. 

In  sec.  7,  it  will  be  seen  that  the  calculation  of  the  theoretical  FLC 
for  plastic  instability  is  based  entirely  on  Y values  less  than  one.  In  order 
to  use  constitutive  equation  parameters  determined  from  tensile  data,  extrap- 
olation to  values  of  f at  c > c is  required.  The  comparison  to  the  ALCOA 
bulge  test  data  supports  such  extrapolation  if  the  procedures  developed  here 
for  modeling  the  tensile  strain  hardening  behavior  are  employed. 

5.2  Load  Relaxation  Tests. 

Lead  relaxation  tests  were  used  to  investigate  the  dependence  of  the  flow 
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stress  o and  the  strain  rate  sensitivity 

v=  (Slno/Jln  i)  (5-7) 

l 

on  strain  rate  and  temperature . For  the  load  relaxation  tests,  the  specimen 
was  loaded  to  a plastic  strain  of  approximately  l'i  at  a strain  rate  of 
10  ^ec  During  loading  and  initial  relaxation,  load  readings  were  taken 
at  a rate  of  appiKiximately  15  per  6econd.  The  data  acquisition  system  was 
programmed  to  reduce  this  rate  incrementally  during  the  tests  which  lasted 
from  2 to  5 hours.  The  load-time  data  was  processed  to  calculate  true 
stress  0 and  strain  rate  > at  a subset  of  load  intervals  separated  by  approx- 
imately equal  load  increments.  The  results  are  presented  as  log  vs.  log 
plots  termed  hardness  curves.  Since  the  strain  accumulated  during  a load  relax- 
ation test  is  very  small,  the  hardness  curve  represents  a specimen  at  essen- 
tially constant  structure,  and  the  slope  of  the  hardness  curve  is  the  strain 
rate  sensitivity  v. 

The  hardness  curve  for  a 2024-0  aluminum  sample  obtained  at  25°f  is 
shown  in  Fig.  14.  The  slope  of  the  hardness  curve  is  essentially  independ- 
ent  of  strain  rate  and  gives  a strain  rate  sensitivity  v = 3\10  This  small 
value  of  v is  in  general  agreement  with  the  lack  of  strain  rate  dependence  of 
the  flow  stress  determined  for  the  constant  extension  rate  tests  at  25"C.  For 

example,  if  a factor  *v  had  been  included  in  Lq.  (5-7),  its  variation  over  the 

-5  -1  -1 

strain  rate  range  of  10  to  10  sec  would  bo  only  2 percent. 

Hardness  curves  determined  on  another  2024-0  specimen  a!  4 different 
temperatures  are  shown  in  Fig.  15.  These  tests  were  all  run  cm  the  same 
specimen  at  decreasing  temperatures . Approximately  1 to  2 percent  plastic 
strain  was  accumulated  for  each  test  prior  to  the  relaxation.  Hence,  each 
hardness  curve  may  represent  a different  hardness  state  and  the  shapes, 
rather  than  the  relative  stress  levels,  provide  the  usotul  informa  Mon.  It 
is  seen  that,  in  general,  the  strain  rate  sensitivity  v increases  with  increas- 
ing temperature  and  decreasing  strain  rate.  At  1P3'Y  and  '25'Y,  the  hardness 
curves  have  an  inflection  point,  and  the  strain  rate  sensitivity  decreases  w’th 
strain  rate  below  a strain  rate  of  approximately  111  ( sec  \ Those  temperatures 
represent  approximately  0.40  and  0.54  ot  the  Melting  point  oi  t ho  alloy,  respec- 
tively, which  is  a temperature  range  where  one  experts  grain  boundary  sliding 
to  contribute  to  the  deformation.  Previous  studies  !2>u,37]  have  shown  that 
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FIGURE  14.  STRESS  RELAXATION  HARDNESS  CURVE  LOG  a vs.  LOG  i 
FOR  2024-0  ALUMINUM  AT  25°C 
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grain  boundary  sliding  can  lead  tc  a hardness  curve  with  an  inflection 
point  at  Intermediate  strain  rates,  and  it  is  thought  that  grain  boundary 
sliding  determines  the  dependence  of  the  strain  rate  sensitivity  on  strain 
rate  at  163°C  and  225°C  for  2024-0  aluminum. 

Since  the  hardness  curve  at  163°C  is  of  interest  in  regard  to  the 
interpretation  of  the  strain  hardening  data  at  this  temperature,  it  is 
shown  expanded  in  Fig.  16.  The  maximum  strain  rate  sensivity  is 
v . 6x1 0 at  r.  ; 10"5  sec~^ , and  v may  decrease  by  up  to  a factor  of  10 
as  the  strain  rate  increases  to  10  ^ sec  V This  is  reminiscent  of  the 
strain  rate  dependence  of  the  flow  stress  and  strain  hardening  coefficient  >, 
but  a detailed  model  is  required  to  make  the  comparison  quantitative.  In 
summary,  the  load  relaxation  tests  have  provided  interesting  information  on 
the  strain  rate  dependence  but  further  development  of  constitutive  equations 
is  required. 


43 


LOG  5TRE55  C MPR > 


6.  RESULTS  FOR  Ti-6Ae-4V 


b - 1 Const ant  Extension  Rate  Tests 

Stress-strain  curves  and  tensile  test  parameters  - sheet.  Constant 
extension  rate  tests  were  performed  on  the  T i - 6 A w - 4 V test  specimens  cut 
from  the  RMI  sheet  described  in  sec.  4.1  at  25°C  and  strain  rates  of  10 
to  10  ^ sec  \ In  addition,  some  tensile  test  results  were  available  from 
work  performed  at  the  University  of  Kentucky  (UK)  [4]  at  25°C,  538°C 
(1000°F),  and  6 7 7 °C  (1250°F).  The  1250°F  temperature  is  within  the  warm 
forming  range  for  Ti-6Ac-4V.  However,  some  preforming  is  done  at  room 
temperature,  and  the  25°C  tensile  data  is  relevant  for  this  application. 

The  25°C  tensile  tests  were  performed  on  0°(L)  specimens  at  strain 

rates  of  10  ^ to  10  ^ sec  ^ and  on  90°(T)  and  45°(F)  specimens  at  strain 
-4  -2  -1 

rates  of  10  to  10  sec  . A typical  Instron  chart  record  obtained  on 

-3  -1 

an  L specimen  at  10  sec  is  shown  in  Fig.  17.  The  load-elongation  curve 
shows  very  little  strain  hardening.  There  is  a mild  yield  point  effect  follow- 
ed by  a slight  rise  to  a maximum  locd.  In  general , the  strain  hardening  was 

-2  -1 

largest  at  the  lowest  strain  rates.  At  10  sec:  and  above,  the  load  tended 
to  fall  off  continuously  after  yielding.  The  thickness  strain  transducer  was 
used  for  these  tests,  and  il  is  seen  that  the  transducer  output  is  smooth  and 
nearly  1 i near. 

The  load-elongation  curves  have  been  used  to  compute  standard  tensile  test 

parameters  vs.  strain  rat".  These  are  presented  in  Table  9.  1 he  engineering 

stress  at  the  yield  point  o is  recorded  instead  of  a 0.2  percent  yield  stress 

Both  the  yield  point  stress  u and  the  ultimate  tensile  stress  cllTr  increase 

u UTS 

steadily  with  strain  rate,  Goth  the  uniform  strain  t and  the  total  elon- 
17  umf 

gation  ■.£  decrease  with  increasing  strain  rate  in  a manner  such  that  the 
post-uniform  elongation  is  nearly  constant.  The  stress  values  are  nearly 
independent  of  orientation  with  those  for  the  T specimens  being  about  2 per- 
cent loss  than  these  for  the  L and  T specimens.  The  values  are  in  excel- 

lent agreement  with  average  values  determined  by  Batelle[4]  (l : 1014  MPa 
T:  1005  MPa)  at  a strain  rate  of  1.67  x 10  ^ sec  . 

The  variation  in  tensile  properties  with  strain  rate  is  also  illustrated 
in  fig.  18  by  the  true  stress-true  strain  curves  (up  to  . -f)  for  the  four 

# o i unit 

L specimens  tested  at  strain  rates  i - 10  ' sec  . It  is  seen  that  the  stress 
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TABLE  9.  TENSILE  TEST  PARAMETERS  FOR 
Ti-6Ae-4V  SHEET  AT  25°C 


" ■ 11  1 ■ 1 J ■*■ 

Specimen 

Rate  (sec~^) 

ou(MPa) 

oUTS(MPa) 

cunif 

,:t0t(5;) 

24L 

10'1 

1064.5 

-- 

-- 

8.2 

25L 

10'2 

1035.6 

1029.1 

0.025 

10.1 

26L 

10'3 

998.1 

1000.0 

0.053 

12.4 

271. 

10‘4 

965.6 

985.7 

0.069 

12.7 

28L 

10'5 

938.1 

965.5 

0.070 

13.4 

17T 

10'2 

1033.0 

-- 

-- 

6.9 

16T 

10'3 

999.0 

996.4 

0.040 

10.3 

18T 

lo'4 

972.5 

983.3 

0.057 

11.8 

1 OF 

10“2 

1018.7 

-- 

-- 

9.5 

17F 

10-3 

387.6 

978.1 

0.025 

13.2 

19F 

10'4 

956.8 

960.0 

0.049 

13.1 

strain  curves  are  nearly  parallel. 

Tensile  data  on  specimens  cut  from  the  same  T 1 -6Ac -4 V sheet  has  also  been 

obtained  at  UK  In  a parallel  study.  Their  true  stress-true  strain  curves  vs. 

temperature  and  strain  rate  are  shown  in  Fig.  19  and  Fig.  20  taken  from 

reference  [ 4 ].  In  general,  the  flow  stresses  decrease  and  total  elongations 

Increase  with  Increasing  temperature  and  decreasing  strain  rate.  The  room 

-2  0 -1 

temperature  stress-strain  curves  at  10  and  10  sec  are  considered  anomalous. 

Stress-strain  curves  and  tensile  test  parameters  - bar.  Since  the 
T1-6AE-4V  sheet  work  hardens  very  little,  some  tensile  tests  on  Ti-6Af-4V 
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bar  performed  for  a previous  study  will  be  reexamined  here.  These  tests  were 
run  on  button-head  specimens  at  a strain  rate  of  10  c sec  and  at  temperatures 
of  25°,  200°,  350°,  and  500°C.  An  Instron  chart  record  obtained  at  25°C  show- 
ing load  vs.  elongation  is  shown  in  Fig.  21  and  the  true  stress-true  strain 
curves  (up  to  i-  for  all  four  temperatures  are  shown  in  Fig.  22.  It  is 

seen  that  the  strain  hardening  is  significantly  greater  for  the  bar  specimens 
than  for  the  sheet  specimens.  In  fact,  the  load-elongation  curves  for  the 
sheet  specimens  look  very  much  like  those  for  bar  specimens  that  have  been 
subjected  to  a 1 or  2 percent  tensile  prestrain  [29  ].  This  is  consistent 
with  the  observation  that  the  sheet  material  has  not  been  fully  annealed. 

The  tensile  test  parameters  for  the  bar  specimens  vs.  temperature  are 
shown  in  Table  10.  At  25‘',C,  the  yield  stress  c is  smaller  and  the  ultimate 
stress  oyyr  is  larger  than  for  the  sheet  material. 

TABLE  10.  TENSILE  TEST  PARAMETERS  FOR 
T i -6A c - 4 V BAR  AT  10'2  sec'1  . 


H"—  ^ r L-  ' " ‘ ' 

Temperature 

ay(MPa) 

Ours (MPa) 

‘'uni  f 

r-f 

O 

ft 

>■> 

25°C 

932.2 

1056.1 

0.045 

15.1 

200°C 

639.8 

S61.3 

0.048 

14.9 

350°C 

626.2 

789.2 

0.054 

14. a 

500"C 

5bb . 0 

674.5 

0.054 

17.9 

Plastic  Anisotropy.  In  the  parallel  study  on  the  1 i -6A c -4V  sheet  at  UK 
[ 4 ],  extensive  work  was  done  on  measurement  of  r-values  vs.  orientation, 
strain,  strain  rate,  and  temperature.  The  r-values  were  determined  by  incas- 
ing changes  in  dimension  of  0.05  inch  square  grids  that  had  been  etched  on 
the  flat  surface  of  the  tensile  specimen  prior  to  testing.  The  UK  results 
at  room  temperature  and  3 percent  strain  are  shown  as  a function  of  orienta- 
tion and  strain  rate  in  Table  11.  The  r-values  range  from  0.49  to  0.30. 

The  scatter  is  quite  large,  and  no  definite  patterns  regarding  dependence 
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FIGURE  21.  LOAD  vs.  ELONGATION  FOR  Ti-6Ae-4V  BAR  AT  25°C 
AND  10  ? sec  1 
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on  strain  rate  or  orientation  emerge.  However,  it  was  shown  that  the  r-values 
Increased  substantially  with  strain,  Increasing  from  '0.3  at  yield  to  '1.0 
at  a true  strain  of  0.1. 


TABLE  11.  R-VALUES  FOR  Ti-6Ae-4V  SHEET 

AT  ROOM  TEMPERATURE.  FROM  [ 4 ]. 


Orientation 

at  10^  (sec-^ ) 

at  10  ^ (sec  ^ ) 

at  10'4  (sec 

ro 

0.51-0.55 

0.58-0.69 

0.60-0.64 

r90 

0.51-0.52 

0.52-0.54 

0.50-0.51 

r45 

0.59-0.80 

0.70-0.73 

0.49-0.58 

Measurements  in  our  laboratory  based  upon  the  transverse  strain  sensor 
output  also  showed  significant  scatter.  Using  Eq.  4-2,  measured  r-values 
ranged  from  0.33  to  0.74  with  no  definite  dependence  on  strain  rate  or  orien- 
tation. It  is  interesting  to  note  (Fig.  17)  that  the  slope  of  the  trans- 
verse strain  sensor  output  decreases  over  the  first  '1  percent  plastic  strain 
and  then  is  nearly  constant  through  the  maximum  load  point.  Eq.  (4-2)  can  be 
written 

r=(^)_17-l  (6-1 ) 

where  dt  is  the  decrease  in  thickness  and  dc  the  increase  in  guage  length. 

The  derivative  dt/dv  is  the  slope  observed  in  Fig.  17.  Eq.  (6-1)  together 
with  the  observed  transverse  strain  sensor  output  imply  that  the  r-value 
should  increase  over  the  first  '1  percent  strain  and  then  decrease  as  the 
ratio  ( t / v ) decreases.  This  is  contrary  to  the  observations  on  the  gridded 
specimens,  and  the  discrepancy  is  not  understood. 

As  a final  comment,  it  is  noted  that  the  yield  point  stress  o and 
ultimate  stress  °yy<jare  nearly  independent  of  orientation  at  fixed  strain 
rate  (Table  9 ),  with  the  stresses  for  the  F specimens  slightly  less 
(1.5  - 2.5  percent)  than  those  for  the  L and  T specimens.  If  the  Hill 
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theory  for  anisotropy  holds,  this  would  imply  that  r^  would  he  i»  'qht.1y 
greater  than  rp  and  rgQ.  This  is  suggested  but  not  clearly  demons t ted 
by  the  measured  revalues. 

In  suninary,  the  plastic  anisotropy  of  the  Ti-GA.-4V  sheet  is  charac- 
terized by  approximate  isotropy  within  the  plane  of  the  sheet  (normal 
anisotropy)  with  an  average  r-value  of  approximately  0.6. 

Strain  hardening  behavior  - bar.  Since  the  Ti-6Ac-4V  sheet  material 
used  in  this  investigation  strain  hardens  very  little,  the  strain  hardening 
behavior  of  the  annealed  bar  will  be  examined  first.  The  method  of  analysis 
will  be  the  same  as  that  used  for  the  2024-0  aluminum  in  sec.  5. 


A plot  of  log  x vs.  log  e for  the  25^0  specimen  of  Table  10  is  shown 
in  Fig.  23.  The  shape  of  the  curve  is  clearly  of  t he  same  form  as  those 
obtained  for  2024-0  aluminum.  Again,  log  x decreases  linearly  with  log  o 
at  large  stress.  Reference  to  the  stress-strain  curve  indicates  that  this 
linear  portion  begins  at  approximately  2 percent  strain.  The  straight  line 
fit  to  this  portion  of  the  log  > vs.  log  o data  is  also  shown  in  Fig.  23. 
Using  the  parameter  for  this  line  and  the  stress-strain  data,  the  Swift 
equation  for  the  large  strain  portion  of  the  tensile  stress-strain  curve 
is 

o (MPa)  = 1239  (-0.008  a t )°;03b  (6.2) 


It  is  emphasized  that  this  equation  only  holds  at  <-5"C  and  a strain  rate  or 

.9.1 

10  ^ sec  . The  strain  rate  dependence  at  25°C  has  not  been  investigated. 
The  stress-strain  curve  and  Swift  equation  fit  are  shown  in  Tig.  24.  The 
predicted  uniform  strain  n-co  = 0.043  agrees  closely  wit.li  the  measured  value 


of 


u n i f 


= 0.045. 


The  portion  of  the  log  ^ vs.  log  o plot  corresponding  to  strains  from 
approximately  0.005  to  0.020  (log  'i  = 0.5  to  1.0)  can  also  he  fit  by  a 
straight  line.  The  strain  hardening  exponent  for  this  portion  is  n = 0.06. 
For  strains  less  than  0.005,  the  log  a vs.  loo  o data  rises  above  this 
straight  line  due  to  a transient  strain  contribution. 

The  strain  hardening  behavior  of  the  Ti-6A>-4V  bar  was  also  investi- 
gated for  the  elevated  temperature  runs  in  terms  ot  log  vs.  log  . 


Although  these  plots  had  generally  the  same  shape  as  that  for  25°C,  it 
was  more  difficult  to  identify  a high  stress  linear  portion.  This  is 
illustrated  by  the  plot  for  the  500JC  tensile  test  in  Fig.  2b.  Since 
the  identification  of  a straight  line  segment  is  not  clear,  the  modeling 
procedure  for  the  high  temperature  data  will  not  be  reported. 


It  is  significant  that  the  strain  hardening  behavior  characterized 
by  a linear  dependence  of  log  -y  on  log  e has  been  observed  for  both 
2024-0  aluminum  and  Ti-6Ac.-4V,  alloy  systems  with  significantly  different 
structure  and  properties. 


Strain  hardening  behavior  - sheet.  It  was  possible  to  analyze  the  strain 
hardening  behavior  for  the  Ti-6A(-4V  sheet  material  in  terms  of  log  vs. 
log  o for  those  tensile  tests  that  exhibited  a maximum  load.  These  results 
will  be  sunniarized  for  the  four  L specimens  tested  af  strain  rates  of  10  ^ 
to  10'"  sec'1  . 

-4  -1 

The  log  a vs.  log  c plot  for  the  test  run  at  10  sec  is  shown  in 
Fig.  26.  It  is  essentially  a straiont  line  over  the  entire  range  of 
stresses.  The  points  begin  at  a strain  c = 0.015  chosen  to  avoid  the  region 
of  minimum  stress  which  follow  the  yield  point.  Straight  lines  were  fit  to 

the  log  i vs.  log  c data  corresponding  to  strains  greater  than  2 percent 

-2  -1 

(1.5  percent  for  the  10  ‘ sec  test).  The  stress-strain  curve  and  the 
Swift  equation  fit  for  the  10  ^ sec  1 test  is  shown  in  Fig.  27.  The  fit 


is  excellent  over  the  entire  range  of  strain. 

The  Swift  equation  parameters  for  the  four  L specimens  are  summarized 
in  Table  12.  The  n values  are  nearly  independent,  of  strain  rate  with  an 
average  value  n=0.106.  This  is  higher  than  either  the  low  strain  (0.06) 
or  high  strain  (0.035)  values  for  the  annealed  bar.  The  ■.  values  increase 
with  strain  rate  sc  as  to  account  for  the  decrease  in  uniform  strain.  Again, 
the  agreement  between  calculated  and  measured  uni  form  strains  is  very  good. 
The  k values  generally  increase  with  strain  rate. 

In  summary,  the  strain  hardening  behavior  of  the  Ti-6A.-4V  sheet  can 
be  described  by  a single  Swift  equation  over  the  entire  range  of  strain 
past  the  initial  yield  point.  Although  an  average  n value  could  be  used, 
the  i and  k values  depend  upon  strain  rate. 
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TABLE  12.  SWIFT  EQUATION  PARAMETERS  FOR 
T1-6AL-4V  SHEET  AT  25°C. 


Specimen 

Rate 
(sec-1 ) 

n 

t: 

O 

k(MPa) 

n-Lo 

runi  f 

25L 

10-2 

0.093 

0.073 

1310 

0.020 

0.025 

26L 

10~3 

0.121 

0.063 

1367 

0.058 

0.053 

27L 

10-4 

0.102 

0.033 

1334 

0.069 

0.069 

20L 

O 

1 

LT 

0.106 

0.031 

132? 

0.075 

0.070 

6 . 2 Load  Relaxation  Tests. 

Variation  witli  tempo nature  - bar.  Load  relaxation  tests  have  been  per- 
formed to  determine  the  dependence  of  the  flow  stress  and  strain  rate  sensi- 
tivity on  strain  rate  and  temperature  for  the  annealed  bar  material.  Six 
load  relaxation  tests  were  performed  on  a simile  specimen  at  decreasing 
temperatures  of  500°,  425° , 300°,  275°,  200",  and  25"C.  The  hardness  curves 
are  shown  in  Fig.  28  and  the  four  runs  at  the  intermediate  temperatures  are 
shown  on  a more  expanded  scale  in  Fig.  29.  These  hardness  curves  illustrate 
the  change  from  low  temperature  to  high  temperature  grain  matrix  deformation 
behavior  suggested  by  Hart  [16].  For  the  strain  rates  sampled,  the  transi- 
tion occurs  at  35P'T..  At  temperatures  below  350"f,  the  strain  rate  sensi- 
tivity is  small  and  increases  with  increasing  strain  rate  at  each  temperature. 
A typical  value  would  be  v , 0.01  at  ?6"C  and  10  4 sec  \ At  temperatures 
above  380° C,  the  strain  rate  sensitivity  is  lame  at  low  strain  rates  and 

-i,  _i 

decreases  with  increasing  strain  rat'1.  At  SOOT  and  It)  sec  , it  is 

v 0. 11).  At  350'C.  the  strain  rate  sensitivity  passes  through  a minimum 
-3  -S  -1 

value  of  v . 2x10  " at  a strain  rate  of  10  ' sec  . 11  is  expected  that 

each  of  the  hardness  curves  in  fig.  28  could  hr  til  hy  Ig.  (?-(>)  for 

constant  values  of  \ and  M.  liie  parameters  .-■*  am!  . ° would  increase  as 

• ^ 

strain  is  accumulated  from  run  to  run  and  i would  carry  the  major  tempera- 
ture dependence. 
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Variation  with  strain  - bar.  Two  load  relaxation  tests  were  performed 
on  another  specimen,  both  at  25°C,  at  true  strains  of  0.015  and  0.051.  The 
hardness  curves  are  shown  in  Fig.  30. 

At  25°C,  Hart's  equation  Eq.  ( 2-6)  relating  flow  stress  to  strain  rate 
can  be  approximated  by 

c = o*  + aQ  e (6-2) 

This  equation  was  fit  to  the  two  hardness  curves  in  Fig.  30.  A value  M=30 
was  selected  by  trial  and  error,  and  a*  and  oQ  were  then  determined  by  linear 
regression.  These  are  shown  in  Table  13.  The  value  of  M should  be  con- 
sidered approximate.  Its  choice  affects  the  absolute  values  of  a*  and  oQ 
but  not  the  change  from  test  to  test.  It  is  seen  that  o*  increases  with 

strain  but  n is  almost  constant, 
c 

TABLE  13.  HARDNESS  CURVE  PARAMETERS 
FOR  Ti -6At-4V  AT  25°C. 


Specimen 

Test  No. 

M 

o*(MPa) 

o0(MPa-sec1/M) 

# 2 (Bar) 

1 

30 

683.6 

418.3 

# 2 (Bar) 

2 

30 

740.5 

428.6 

#39L  (Sheet) 

1 

30 

611.1 

490.9 

Comparison  of  load  relaxation  and  constant  extension  rate  tests  - sheet. 

A single  load  relaxation  test  was  performed  on  a Ti-6Ai-4V  sheet  specimen  at 

-3  -1 

25°C.  The  specimen  was  loaded  at  an  initial  strain  rate  of  4x10  sec  and 
the  relaxation  began  at  a strain  of  0.016  and  a stress  level  nearly  equal  to 
the  yield  point  stress  o . Since  the  sheet  material  strain  hardens  so  little, 
a repeat  test  on  the  same  specimen  is  riot  of  interest. 

The  hardness  curve  for  this  test  is  shown  in  Fig.  31  along  with  the 
fit  of  Eg.  (6-2)  for  M=30.  The  a*  and  o values  are  included  in  Table  13. 
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It  is  seen  that  o*  is  smaller  and  c*  larger  than  chose  for  either  test  on 
the  bar  specimen.  Thus  it  is  u rather  than  o*  which  correlates  with 
relatively  high  yield  point  stress  of  the  sheet  material. 

The  dependence  of  flow  stress  on  strain  rate  from  the  load  relaxation 
tests  can  also  be  compared  to  the  strain  rate  dependence  determined  from 
the  constant  extension  rate  tests.  For  this  comparison,  the  yield  point 
stresses  o from  Table  9 were  converted  to  true  stresses  and  plotted  vs. 
initial  strain  rate  in  Fig.  31.  It  is  seen  that  the  load  relaxation 
results,  including  the  extrapolation  to  10’  sec'1,  agree  extremely  well 
with  the  yield  point  stresses  from  the  constant  extension  rate  tests  at 
all  five  initial  strain  rates.  It  is  believed  that  this  close  agreement 
is  due  to  the  fact  that  the  strain  hardening  rate  is  small  and  nearly 
independent  of  strain  rate. 


For  the  Ti-6Ac-4V  sheet  material,  it  appears  that  an  approximate 
strain  and  strain  rate  dependent  constitutive  equation  can  he  obtained 
by  combining  Eg.  (6-2)  with  the  Swift  equation.  A possible  form  would  be 


a = o*  (1  + 

where  o*  and  a are  taken  from  Table 
o 

would  only  be  valid  for  strains  from 


1 


0.106.  • 1/M 

+ ci  e 


13,  e from  Table  12. 
o 

0 to  (0.106  - c ). 


(6-3) 

The  equation 
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7.  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  FORMING  UNITS 


For  2024-0  Aluminum  at  25*C.  Both  in-plane  and  out-of-plane  FLD's 
have  been  obtained  for  the  2024-0  aluminum  at  room  tempera  Lure  by  ALCOA 
[4].  The  out-of-plane  FLD  was  obtained  using  a modification  of  Hecker's 
[19]  method,  and  the  in-plane  FLD  was  based  on  both  standard  parallel - 
sided  and  noticed  tensile  specimens.  For  the  notched  tensile  specimens, 
a special  square  grid  of  100  lines/inch  was  used  to  determine  the  strains. 
The  FLD's  are  shown  in  Fig.  32  for  two  different  lots  of  2024-0  aluminum, 
the  standard  Reynolds  lot  and  an  ALCOA  lot.  The  ALCOA  lot  is  considerably 
more  formable  than  the  Reynolds  lot,  and  for  both  lots,  the  in-plane  forming 
limit  is  well  below  the  out-of-plane  forming  limit.  Possible  reasons  for 
this  latter  discrepancy  have  been  discussed  in  Section  3. 

The  Hill  FLC  was  calculated  for  the  Reynolds  lot  2024-0  using  the 
average  Swift  equation  parameters  from  Table  7 and  an  average  r - C.73. 

The  FLC  is  shown  in  Fig.  33  along  with  the  in-plane  and  out-of-plane 
FLD's  for  the  Reynolds  lot  (converted  to  true  strains)  from  Fig.  32.  The 
agreement  between  the  FLC  and  in-plane  FLD  is  very  good.  This  is  due  in 
part  to  the  fact  that  the  strain  rate  sensitivity  of  2024-0  aluminum  at 
room  temperature  is  very  low  (Section  5.2)  and  its  neglect  by  using  the 
Hill  theory  is  a good  approximation.  It  is  not  surprising  that  the  FLC 
agrees  with  the  in-plane  rather  than  out-of-plane  FLD  since  the  FLC  was 
determined  from  parameters  measured  under  plane  stress  conditions.  The 
calculated  FLC  for  2024-0  aluminum  at  25<’C  is  discussed  in  more  detail  by 
Nagpa 1 et  al . [38] . 

For  2024-0  Aluminum  at  163°C.  From  the  Swift  equation  parameters  in 
Table  3 it  can  be  determined  that  the  calculated  Hill  FLC  for  163rC  would 
Mil  below  that  for  25<1C.  The  plane  strain  intercept  fur  163 'C  would  be 
0.122  (assuming  r = 0.73)  vs.  0.155  at  25°C.  This  is  clearly  not  realistic 
as  we  expect  considerably  more  formability  at  163°C  than  at  25°C,  so  the 
FLC  is  not  presented.  It  does  illustrate  that  the  increase  in  formability 
with  temperature  would  be  a strain-rate  effect. 
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TRUE  PLR5TIC  STRAIN  E2 

FIGURE  33.  THEORETICAL  HILL  FLC  FOR  2034-0  ALUMINUM  AT  25l'C 
WITH  COMPARISON  TO  EXPERIMENTAL  FLD’s 


71 


TRUE  PI.R5TIC  STRAIN  E! 


For  Ti-6Av-4V.  Both  in-plane  and  out-of-plane  FLD's  were  determined 
for  Ti-6Ac-4V  at  24°C,  538°C  (1000°F),  and  677°C  (1250°F)  by  the  University 
of  Kentucky  [4].  The  results  are  shown  in  Fig.  34.  For  these  tests,  the 
in-plane  and  out-of-plane  FLD's  agreed  very  closely.  It  is  seen  that  the 
minimum  limit  strain  at  24°C  is  about  0.13  and  the  fomiability  increases 
with  temperature. 

If  a Hill  FLC  were  calculated  based  on  the  Swift  equation  parameters 

in  Table  12,  the  overall  level  would  depend  upon  strain  rate.  The  effective 

strain  rate  for  the  FLD  can  be  estimated  from  the  crosshead  speed  to  be 
-2  -1  -2  -1 

> 10  sec  . Using  the  10  sec  parameters  would  give  a plane  strain 
intercept  of  approximately  0.02.  Again,  this  indicates  that  the  formability 
of  the  Ti-6Av-4V  is  mainly  strain  rate  dependent,  and  the  Hill  FLC  will  not 
predict  the  forming  limits. 


Ti  - 6 Al  - 4V 

Crosshead  Speed:  25in/min(l.06cm/s) 


8.  CONCLUSIONS 


1.  The  results  of  this  program  suggest  that  sheet  metal  formability 
can  be  predicted  from  theoretical  forming  limit  curves  based  on  constitutive 
equation  parameters  determined  from  uniaxial  tension  tests.  Success  of 
this  approach  requires  careful  modeling  of  strain  hardening  behavior  and 
the  consideration  of  strain  rate  sensitivity.  Further  work  to  incorporate 
strain  rate  effects  into  calculated  FLC's  is  required. 

2.  A new  method  of  investigating  strain  hardening  behavior  based 
upon  the  dependence  y - y (a)  is  proposed.  The  method  assures  that  the 
slope  of  the  stress-strain  curve,  which  directly  effects  formability, 
will  be  accurately  modelled. 

3.  Analysis  of  the  stress  dependence  of  the  strain  hardening 
coefficient  y (a)  has  shown  that  the  strain  hardening  behavior  of  2024-0 
aluminum  at  25°C  and  163°C  and  of  Ti-6A.t-4V  at  25°C  can  be  modelled  by  the 
Swift  equation  for  large  strains  (e  * 0.05  for  2024-0  At  and  e s 0.02  for 
Ti-6At-4V).  Excellent  agreement  is  obtained  between  measured  and  calculated 
uniform  strains.  For  2024-0  At  at  25°C,  good  agreement  is  obtained  between 
the  extrapolated  Swift  equation  and  hydraulic  bulge  test  data  to  strains 

of  "0.4. 

4.  The  strain  rate  dependence  a (e)  from  load  relaxation  tests  at 
25°C  can  be  modeled  by  an  equation  suggested  by  Hart  for  low  homologous 
temperatures.  For  2024-0  At  at  25°C,  the  strain  rate  sensitivity  is  small 
enough  to  be  neglected.  For  Ti-6At-4V  at  25°C,  good  agreement  is  obtained 
for  o (c)  determined  from  both  load  relaxation  tests  and  constant  extension 
rate  flow  stress  vs.  strain  rate. 

5.  The  strain  rate  dependence  of  the  flow  stress  a (e)  at  elevated 
temperatures  (to  225°c  for  2024-0  At;  to  500°C  for  Ti-6At-4V)  shows  the 
transition  from  low  temperature  to  high  temperature  behavior  suggested  by 
Hart.  For  2024-0  At  at  163°C  (325°F),  it  appears  that  the  strain  rate 
dependence  is  effected  by  grain  boundary  sliding. 

6.  The  orientation  dependence  of  the  flow  stress  of  2024-0  aluminum 
sheet  at  25°C  can  be  interpreted  using  Hill's  theory  of  planar  anistropy 


and  measured  r-values.  The  agreement  with  effective  stress-strain 
determined  from  the  hydraulic  bulge  test  is  reasonable. 

7.  For  2024-0  Ao  at  25°C,  the  theoretical  FLC  based  on  Hill's 
theory  of  localized  necking  (strain-rate  independent)  is  in  good  agreement 
with  the  experimental  in-plane  forming  limit  diagram.  However,  both  are 
significantly  lower  than  the  standard  out-of-plane  FLO  determined  by 
punch-stretch  tests. 

8.  For  2024-0  aluminum  at  163'T  and  Ti-6A,-4V  at  2B''C,  the  forming 
limits  calculated  using  Hill's  theory  of  localized  necking  would  be  too 
small  due  to  the  neglect  of  strain  rate  dependence.  Development  of  a 
strain  rate  dependent  model  for  the  calculated  FLC  is  required. 
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APPENDIX  A:  SUMMARY  OF  SOME  ANISOTROPIC  PLASTICITY  RELATIONSHIPS 


Hill  [27]  has  extended  the  Von  Mises  yield  function  to  account  for  plas- 
tic anisotropy.  The  resulting  yield  function  f(o..)  and  effective  stress  o 

' J 

are  defined  by 

2 f (o . . )=  F(o  - a ) 1 + G(a  - a )2  + H(o  - o )2 
ij  ' ' y z 7 v z x7  ' x y 7 

^L,y-  .2N,x; 

= (2/3) (F  . G * H)  (A1 ) 


where  F,  G,  H,  L,  M,  N are  the  anisotropy  parameters.  For  the  case  of  planar 
stress  with  = t^2  = 0,  Eq.  (A1)  reduces  to 

2f(o  (G  + H)  ox2  - 2 H + (F  + H)Oy  + 2 N = 2k^2  (A2) 

where  = (F  + G + H)/3.  For  application  to  the  deformation  of  a thin  sheet, 
the  \ direction  is  taken  parallel  to  the  rolling  direction. 

Lee  and  Kobayashi  [2  ] have  presented  the  corresponding  expressions  for 
an  applied  stress  system  defined  with  respect  to  orthogonal  axes  (1,2)  which 
are  rotated  by  an  angle  0 with  respect  to  the  axes  (x,y).  For  principal 
stresses  (il2  = 0),  the  result  is 


?f(o,i)-  Aoj 2 - 2 BoiO:  + Co2?  = 2KAo2  (A3) 

where  A * F$in20  + Gcos20  + H - fsin20  cos2o 
B = H - ic.sin20  cos2o 
C = G sin:0  Fcos ?0  + H - »Ssin2o  cos'O 
6 = F + G + 4H  - 2N. 


The  normal  plastic  strain 
cording  to 


dr  i 


dr  3f 
2kflo 


increments  are  derived  from  the  yield  function  ac- 
* (A4) 


dr; 


JL_  Af 

2kAn 


(A5) 
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The  normal  strain  increment  de3  can  be  obtained  from  the  condition  for  con- 
stant volume 


de3  + de2  + de3  = 0.  (/\g) 

It  is  important  to  note  that  the  shear  strain  increment  dyi2  is  not  necessar- 
ily zero  even  for  a stress  system  with  n2  = 0.  In  fact,  the  principal  axes 
of  stress  and  strain  will  only  coincide  if  one  of  the  principal  stress  axes 
lies  along  the  rolling  direction  [2]. 

In  the  following  we  denote  the  ratio  of  principal  stresses  a2/oi  = a, 
and  the  corresponding  ratio  of  plastic  strain  increments  de2/dei  = 3.  For 
the  case  of  proportional  straining,  6 is  a constant,  and  B = e2/d.  Eqs. 

(A4)  and  (A5)  can  be  used  to  relate  a and  6, 

B = (Cot  - B)/(A  - Bet),  (A7) 


and  inverting, 

a = (AB  + B)/(B3  + C).  (A8) 

We  also  denote  the  ratio  of  plastic  strain  increments  de2/de3  = r.  Eqs.  (A6) 
and  (A7)  can  be  used  to  relate  r to  a and  8, 

r = -B/O+B)  = (B-Ca)/[A-B  + (C-B)ct].  (A9) 


Thus  the  r-value  measured  in  a tensile  test  (a  = 0)  for  a specimen  cut  from 
a sheet  at  an  angle  0 to  the  rolling  direction  is  given  by 


rG(ct=0)  = B/(A-B). 


(A10) 


For  the  usual  case  of  specimens  cut  at  0°,  90°,  and  45°  to  the  rolling  direc 
tion,  Eq.  (A10)  gives  the  relations  between  the  measured  r-values  and  the 
anisotropy  parameters 


r0 


r4s 


N 


1 

7 ’ 


(All) 


or,  inverting, 

1=1  1-1  N n +1 ' 

H “ r0  ’ H " r90  ’ H " r0  r9o 

- 

It  is  now  possible  using  Eqs,  (A3)  and  (A12)  to  write  down  the  relations 
between  the  effective  stress  0 and  the  measured  stress  for  certain  simple 
test  configurations  in  terms  of  the  measured  r-values.  For  a uniaxial  tensile 


4 s 


1 

7 


(A12) 


1 
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«st  'o,  • or  - 0)  with  a specimen  cat  at  an  anple  « ■ w 45 

to  the  rolling  direction  we  have 


0 - 0°  : 


0 = 90°: 


0 = 45°: 


-3Lr.^rflr„L'l  1/2 , 
2TF^r„+rcr,D] 


^(r^rnrso) 
2(r0+r90+ror so) 


0 ~ jj4( r p + r 9 o + n or  so) 


1/2 

1 

1/2 


(A13) 

(A14) 

(A15) 


For  a hydraulic  bulge  test,  assuming  balanced  biaxial  stress  at  the  pole 
(01  = 02  ~ Og), 

1 1/2 


B 


3ir^ul _ 

2(ro+nso+ror5o ) 


(A16) 


,.,t  nr  (A5)  the  relation  between  the  effective  strain  in- 
Also,  using  Eqs.  (AA)  or  iA5)’  nhtained  For  the  tensile 

crement  <k  and  the  measured  strain  increment  can  he  ohtained. 

test,  _ _ (Al?) 

dc/dct  = 1 ’ 


(A18) 


and  for  the  bulge  test, 

de/de  3 = (5/bg)  1 • 

nf  the  factors  in  F.qs,  (A13)  to  (A16)  relate  effective  and 
Thus  the  inverse  of  the  tactoi  s 

measured  strain  increments. 

. . ..  cnprial  case  of  normal  anisotropy,  i.e.,  iso 

Finally.  «e  » ' ’ _ 1ndepMdwt  #f  From  W.  <M2) 

tropy  within  the  plane  such  that  r^t 

the  isotropy  conditions  are 

G/H  = F/H  = N/H  - 2 a l/r* 

. A B C in  Fq.  (A3)  become  independent  of  0.  By  choosing 

and  the  parameters  A,  B,  t in  v 

H = r/  (1+r) , tq. (A3)  becomes 

(A20) 


(A19) 


2f(c„>-  o1  - fir  * “ Zb 


h . . (,,ru3(Hr)  Also,  k,  replaces  kA  in  tus.tMl  and  (Ab).  tor 

r:,  Lsl^'Ihe’ratid  of  tie  effective  stress  to  the  aria- 


B2 


y 


iWvIlsJtiJI  IkIn  -nlltilMf J"  ^ '-1  WUHfc Jll!  hii^ jjMitfXt**:--. 


yield  stress  is 


(A21) 


Thus,  o can  be  thought  of  as  the  uniaxial  tension  yield  stress  for  an  iso- 
tropic specimen  (r  = 1).  For  the  hydraulic  bulge  test,  Eq.  (A16)  or  Eq. 
(A20)  gives 


c 


Eqs.  (A21)  and  (A22)  combine  to  give  the  familiar  result 


(A22) 


(A23) 


In  some  treatments,  Eq.  (A20)  is  written  in  the  slightly  different  form 


o i‘ 


2r 

T+r 


OlO; 


0;* 


\‘  • 


(A24) 


In  this  case,  \ is  simply  a tensile  yield  stress  and  differs  from  the  effective 

1 /2 

stress  o defined  here  by  the  factor  (2k j)  1 . 
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11. 

APPENDIX  B.  AUTOMATIC  DATA  ACQUISTION 

SYSTEM  CONTROL  PROGRAMS 

0 : " DVM  RUM": 

3 1 : if  •:  p p i J 1 1 = " 

C-u 

r 2 •*  r 4 ■*  r 3 

i - ———————  c 

Y “ : 9 j b " p r t a 1 1 " 

61 

" c sine  w " : 

2:  dim  If [256]. 

3 2 : " m t-  a p e " : b e e p 

62 

if  r 4 > r 5 5 

SI C 3 0 O ] > D f [3000 

33:  ent  "Store 

9 1. o "csifin" 

] 

on  tape ? “ > J 1 

63 

if  S f [ r 4 . 

3 : dim  C $ 1 5 0 ] > 

34:  if  c a pUMI  " 

r 4 ] = " j “ j 9 1 o " < " 

Jf C32] 

Y " 5 e t o " m stop" 

64 

r 4 + 1 9 r 4 ! 9t  o 

4:  dim  2 [53 

35:  fnp  ” t r a c k 

" c s i n e w " 

5 : " ft  C D t 5 1 S P 0 ’ * 

# " . Y 

65 

"esisen": 

Cl 

36:  enp  "file 

66 

POs(CX.Sf[r2 

6:  0*P4rl3*rl4 

tt" . X 

. 

•2]  ) *r6?  r 2 + 

?:  1 * r 1 O*  r 1 1 

37:  trk  Y5fdf  X 

1 

»r2 

3 : cl  1 ’ i n i t f 

38:  idf  Z [1 ] « 

67 

if  r 6 > O 5 9 1 o 

9 : on i ?. "CS I O" 

Z [2 3 » 2 [3]  . Z [4]  . 

•• 

: s i e x " 

10:  wrt.  " d v m " . 

ZC53 

6S 

p r t "CSI 

"F1R3T2M3H008D0 

3S : if  Z C2] #0 : 

ERROR" 5 dsp  "ERR 

e n t " O v e r w rite 

OR!  " 

lls  fmt  8." 

f i 1 e??“ . Jl 

69 

p r t " r 1 - r 6 

40:  if  Z[2]#0: 

- 

" . r 1 . r 2 . r 3 > 

12:  " m 1 o o p " : s p c 

if  cap ( Jf ) # " Y ” ! 

r 4 > r 5 . r 6 

5 5 w r t 1 6 . S 

9 1 o " m rape " 

70 

" c s i f i n " : 1 - C 

13:  9Sb  "I DENT" 

41 : if  Z [4] <3600 

t 

5 

i r e t- 

14:  s p c 2 5 w r t. 

< beep. ent  “File 

71 

" c s i e x’ " : 

1 6 . 8 5 p r t " S e t t i 

too  small" j J#: 

-!•  -A 

1'  i. 

o m p r 6 

n 9 s : " 

9 1 o " m t a p e " 

73 

? t o " R " 

15:  9sb  "SETTING 

42:  ref  X. If. SI. 

74 

9 1 o " C " 

S" 

Of 

75 

9 t.  o " 0 " 

16:  wrt  16. S 

43:  p r t "file 

i 6 

9 t o " . " 

17:  d s p "Hit 

stored" 

~ * 

i*  »■ 

9 1 o " : " 

continue  w ft  e n 

44:  " m s t o p “ : p r t. 

78 

9 1.  o "space" 

ready  to  eta r t . 

"Done” < spc  55 

79 

9 1.  o " I " 

St  P 

SO 

9 t 0 " S " 

IS:  stp 

45:  9 t o " m 1 o o p “ 

SI 

9 t 0 " P " 

IS:  -2*T 

46:  "clock": r 1 7 

82 

9 t 0 " 0 " 

20:  eir  7? wrt 

X ! 0 •*  r 1 7 

83 

" R " : 

"elk". "POOIEISR 

47:  w r t.  " c lk"« 

84 

9Sb  "NUN” 

"TIOOEIDR" 

85 

X - 1 ■>  r ! 

21:  " p t e s t " : e l r 

48:  fmt.  “P".f3.0 

S6 

r 2 4 r 3 

7 

. "E"»f 1.0. "SR” 

37 

9 1 c " C S I " 

22:  if  r 1 7 # 0 ? 

4 9 : w r t "cl k " . 

88 

L • 

9 1 o "clock" 

r 1 6 . X 

S9 

9sb  "MUM” 

23:  if  r 1 3 < 1 e6 : 

50:  eir  7 

90 

c 1 1 J scnlMI: 

9 1 o "pwait" 

51:  9 1 o " p t e s t " 

) 

24:  <9sb  "Print" 

52:  " C 8 1 0 " : 

91 

X9  r 7 

25:  "pwa  it":  if 

53:  T + 1-»T 

92: 

l+r7/10O9Y 

rl7#0$9to  "cloc 

54:  fmt  f 5 r e d 

93: 

Y 9 1 e 6 2 9 Y 

k" 

"dvm" « fl 

94: 

fts  ( Y ) 9 0 f [ r 

26:  if  D#l?9to 

55:  r7-»c 

16. r 1 O + 3 ] 

" p t e s t " 

56:  " C S I ” : 

95: 

r 1 0 + 4 9 r 1 O 

27:  U64rl3 

57:  if  r2>  = r 5 * 

96: 

9 1 0 "CSI" 

28:  if  P # 1 ? 9 1 o 

if  r 1 < = 0*  9t  o 

97: 

" D " : 

"Ptest  ” 

“csif.n" 

98 : 

fmt  "Channel 

29:  beep 

58:  if  r 2 < = r 4 5 

" . f 2 . O . " ) = 

30:  ent  "Print 

9 1 o "esisen" 

**  9 

f 8. 5 

out  all  data?" » 

59:  if  r 1 > 0 J r 1 - 

99 : 

wrt  O . C * 8 

.'t 

1 ■»  r 1 » r 3 ■»  r 2 • ?to 

1O0:  eto  "CSI" 

CS  I •' 


84 


101 : "."Si  ret 

102*  "5  " : 3t  0 

"CSI" 

103s  ''  space"  s 3 t o 

"CSI" 

104:  "I": 

105:  red  "clk“>U 
106:  if  r 14  = 0 5 - 
141.1?  1 4 r 1 4 
107:  <?sb  "HUM" 
103:  X*le3*r8*r 1 
3 

109:  i n t (1 o ? ( r 3 ) 
1 4X 

110:  r 8 / 1 0 T X 4 r 8 
111:  r 8 * 1 9 0 4 r 8 5 
X- 

112:  f nt  "P"? 
f 3 . 0 j " E " j f 1 . 0 f 


SR" 


113: 

r 8 + r 1 6 5 X 4 r 1 

114: 

10+  rS  + X4Y 5 

1+Y/l 00004 Y 

115: 

Y * 1 e 6 1 4 Y 

116: 

fts  (Y)4D$C 

r 1 0 j 

r 10  + 33 

117: 

fts  (U+1)4D 

i Crl0  + 4> r 10  + 73 

1 IS: 

fts  (T)4D*C 

r 10  + 3? r 10+1 13 

119: 

r 1 0 + 1 2 4 r 1 0 

120: 

0 4 T 

121 : 

i ret 

122: 

"S": 

123: 

fts  (fll+DiC 

r 1 0 

» r 1 0 + 33 

124: 

r 10  + 4 + r 10 

125: 

3 1 o "CSI" 

126: 

" P " : 

127: 

fts  ( le63) 4 

Di  C r 1 0 » r 10  + 33 

128: 

r 10  + 44 r 10 

129: 

3 1 o "CSI" 

130: 

" 0 “ : 

131: 

fts  (2e63)4 

0i  Crl0> r 10  + 3] 
132:  rl0  + 4-tr  lO 
133:  gto  “CSI” 
134:  "HUM" s 
135:  if  SiCr2> 
r23  # " ("»< no 
“nunerr" 

136:  r2+l->r2 
137:  sf  3 IS  val( 
SiCr2>  r43  i+V. 
133:  for  J=r2 
to  r4 


139:  if  Sf[J. 

J ] = “ \ " ; 3 t o " n u n 
ret" 

140:  next  J 
141:  “rumerr": 
142:  prt  " n u « 
error!  " 5 d s p 
"ERROR!  " 

14  3:  p r t " r 1 - 
r6"i  r 1 > r 2 > r 3 » 
r 4 > r 5 > r 6 


144: 

140 

145: 

" n u n r e t " : 

146: 

J+  l->r2 

147: 

ret 

148: 

"SETTINGS": 

149: 

••  "4Si 

150: 

prt 

151 : 

04rl4r44r64 

D4X 

152: 

l+r24r34r5 

153: 

ent  "Enter 

Sett  in 3 Con n o. n d 

- ■ " 
•.  • 

« Si C r 53 

154: 

p r t Si 

155: 

gto  +2 

156: 

"set  1 o o p " : e 

rip 

line  = » 

Si C r 53 

157: 

if  Si Cr5 j 

r53 

= "i“  5 spc  • 

prt. 

"len  Si  = "» 

r 5 J 

3to  "set end" 

153: 

ler.  (Si)  + 

1 4 r5 

159:  ” 5 "4Si Cr5? 

r53 5 r5+l4r5 
160:  gto  "setloo 
p " 

161:  "setend":sp 
c 5 n r t " S i " > S i 5 
s p c 2 5 r e t 
162:  "init": 

163:  if  p0=05 
r e n 7 5 c 1 r 7 
164:  fnt  9>f5 
dev  " c 1 k " > 7 1 0 4 p 
8 > "ct  r"  > 725-tp9 < 
" d v m " > 7 2 2 4 p 1 0 ? 
"pt  r"  > 7154pl  1 
165:  dev  "sen" » 
709+p12» "sent" » 
7104p13> " s c n 2 " > 
7 1 1 4p l 4 j "scn3"> 
712+P15 

166:  dev  “svn"» 
7244P16! 8*p7 


167:  1.1 1,  b p p 7 ! 
r d s ( 7 > p 3 j p 4 > 
p 5 ) 4 p 6 5 r 0 + 2 1 ( p 7 
- 3 ) b i t ( 2 » p 5 ) 4 r 0 
: 1 +p74p7 


168: 

if  p7 < 1 7 5 

3 1 O 

-1 

169: 

wtb  "ptr.9" 

< 2 7 » 6 9 

170: 

ret 

171 : 

“sen  1 " : 

172: 

fcit  f 3 . 0 j 

w r t 

"sen" » Pi 

173: 

ret 

174: 

“print": 

175: 

if  rll+3>=r 

10; 

I4p; ret 

176: 

stf (Di [rl  1 » 

r 1 1+33 )4Y 

177: 

rl l+44rl  1 

173: 

int ( los ( Y ) J 

4 X 

179: 

Y/10TX4Y 

130: 

if  X = 6 3 5 

3 1 O 

“Pf 13" 

181: 

if  r 12=0 5 

ret 

132: 

if  X <60 ; 

f xd 

6 ; p r t Y 4 

10TX5 ret 

133: 

if  X=6 1 5 

3 1 O 

"int" 

184: 

f- 

-*> 

ii 

|T*» 

ro 

3 1 0 

"chnl " 

185: 

ret 

186: 

"int": 

137: 

f xd  0 

133: 

prt  "Clock 

ticks® ">stf (Oil 

r 1 1 

? rl 1+33 ) 

139: 

prt  "#Inter 

rupt s= " ? st  f (DiC 

rl 

+ 4 < r 1 1+73 ) 

190 

rl 1+8  + r 1 1 

191 

Y - 1 4 Y 

192 

10004 Y 4Y 

193 

int (Y)4X 

194 

(Y-X) 4104Y 

195 

X 4 1 0 1 Y 4 Y 

196 

fit  3? SPC  2 

197 

prt  "Int. 

= 

* Y 

193 

SPC 

199 

ret 

260 

"chnl": 

201 

fxd  0 

202 

Y - 1 4 Y 

203 

Y 4 1 0 0 4 Y 

85 


204:  spc  2?prt 
"Channel  #">Y 
205:  ret 

206:  "pfl9": 

207:  prt  "print 
208:  if  Y= 1 ? prt 
“on" 5 l*rl2 
209:  if  Y=2>  prt. 
"of f "50*rl2 


210: 

ret 

211: 

"I  DENT": 

212: 

M *»  ^ 

If? " “*Jf 

5 1 4 x 

213: 

ent. 

"Enter 

run 

I.  Li 

. " f Jf 

214: 

prt 

“Run 

I.D. 

" 

215: 

9 1 O 

+ 4 

216: 

" i d 1 oop  " : 

217: 

••  ••  4 

Jf 

218: 

ent 

" line:  " 

> Jf 

219: 

if 

len ( Jf ) > 

165  dsp 

"line 

too 

1 on 

9 " 5 b e e p 5 

wa  i t 

300O5  9 1 o - 

220: 

p r t 

Jf 

221 : 

if 

lenCJf )> 

0?  i f len  C J*)  <16 
5Jf$<"  "*Jf5dMP 
0 

222:  J$  Cl  >163  4 If. 

CX3 > 1 en (I f ) + 1 *X 
223:  if  lenUfl) 
0 5<?to  "idloop" 
224:  wrt.  16.8 
225:  ret 
226:  " p r t a 1 1 " : 
227:  wrt  16.8? 
prt  "Data  Set : " 
5 spc  3 

22S : 0-^P  5 1 ■»  r 1 1 
229:  if  P#i;i*rl 
2>9sb  "print" 
230:  if  Pttl>?to 
-1 

231:  ret. 
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0: 

“DVM  DRTR": 

1 • 
2: 

din  I i [ 2 5 6 3 > 

Si [3003 > Di C30GO 

3 s 

din  Hi  [503 

4: 

“RCOi  5 ISP0"  + 

Hi 

5: 

din  Ci  [2»  2360 

3 

6: 

din  Ji  [53 > 

2 t53 

? • 

1 • 

din  T [53  » P [ 5 3 

» Q [53 

3: 

spc  3 

9: 

prt  “**##**** 

********" 

10: 

prt  "DM2" 

11: 

prt  “******* 

*********“ 

12: 

spc  2 

13: 

1+R  + B 

14: 

0+rl+r4+r6+r 

7 + r 1 3 

15: 

I*r2*r3 

16: 

0 + C 

17: 

1 +Q [ 1 3 +Q  [23 

IS: 

enp  "Tape 

track  #"»T 

19: 

enp  “Data 

File  # " » F 

20: 

t rk  T 5 ldf  F, 

It 

> Si  j Di 

2 1 : 

1 en ( Si ) + r5 

22: 

enp  "Load 

Channe  1 : " > P [ 1 3 

23: 

enp  “Other 

channe 1 : " > P [23 

24: 

prt  " 

25: 

prt  "File 

ID 

M 

26: 

prt  ” 

•i 

27: 

spc  5 prt 

m 

28: 

prt  It 

29: 

prt  " 

> SPC 

2 

30: 

""+Ci[13+Ci[ 

23 

31: 

0 + T [ 1 3 +T  [23  + 

T [ 

33  *T [43  +T  [53 

32: 

prt  "CSI" 

33:  "CSI": 

34:  if  r2>  = r55 
if  rl<=05gto 
"csif in- 
35:  if  r2<=r45 
gto  "csiscn" 

36:  if  r 1 > 0 * rl- 
1 + rl I r3+r25  gto 
- C S I " 

37:  r2+r4+r3 
38:  “csinew": 

39:  if  r4>  r 5 5 
gto  "csifin" 

48:  if  Si  [ r4  > 
r4]  = " 5 " 5 gto 
41:  r4+ 1 + r4 5 gt  o 
"csinew" 

42:  “csiscn”: 

43:  pos (Hi>  Si [r2 
> r 21 ) +r65  r2  + 
l + r2 

44:  if  r6>05gto 
"csiex" 

45:  prt  “CSI 
ERROR" 5 dsp  "ERR 
OR!  " 


46 

prt 

” r 1 - r 6 

= 

" » r 1 > 

r2?  r3j 

r 

r5»  r6 

47 

gto 

"csifin 

48 

"csiex" : 

49 

dnp 

r6 

50 

gto 

"R" 

51 

gto 

"C" 

52 

gto 

"D" 

53 

gto 

" » " 

54 

gto 

••  f M 

55 

gto 

"space " 

56 

gto 

"I" 

57 

9t  0 

"S" 

53 

gto 

"P“ 

59 

gto 

"0" 

68 

"HUM 

*•  • 

• 

61 

if  Si [ r2 > 

r23#"  ("  5 gto 
"nucierr” 

62:  r2+ 1 * r2 
63:  sfg  145 val(S 
i lr2,  r4]  ) + X 
64:  for  J = r 2 to 
r 4 

65:  if  3iEJ»J3=" 
)“5gto  "nun  ret" 
66:  next  J 
67:  "nucierr": 


68:  prt  "nun 
error! " 5 dsp 
"ERROR! " 

69:  prt  "rl-r6"> 
r 1 * r2»  r3>  r4>  r5> 
r6 

70:  "nunret": 

71:  J+l+r25ret 
72:  "DRTR"  : 

73:  if  R> len (Di) 
5 gt  o “outerr" 
74:  stf(Di[fl,R+ 
33  )+D 
75:  R + 4+R 
76:  if  0=050+25 

9t0  +2 

77:  int ( log (abs ( 
D) ) ) +2 

78:  D/10t2+Y 
79:  ret 
80:  "R": 

81:  gsb  "HUM" 

82:  X-l+rl 
83:  r2+r3 
84:  g t o "CSI” 

85:  "C": 

36:  esb  "HUM" 

37:  0+L 
88:  if  X = P C 1 3 5 
1+L 

89:  if  X = P [23  5 
2 + L 

96:  gsb  "DRTR" 
91:  if  2#62  5 gt  o 
"Cerr" 

92:  Y-l+Y 
93:  Y* 1 00+Y+C 
94:  if  C#X5  et  o 
"Cerr" 

95:  gto  "CSI" 

96:  "Cerr": 

97:  prt  "Channel 
Error" 

98:  prt  “Connand 
= ".X 

99:  prt  "Data 
= ">  Y 

100:  gto  "CSI" 
101:  "I": 

102:  esb  "DRTR" 
103:  if  2#615 
gto  "Ierr" 

104:  Y-l+Y 
105:  Y* 1 000+Y 
106:  int ( Y) +2 
107:  (Y-2)*10+Y 
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IDS: 

Z> lOt Y+Y 

109« 

Y/U6  + Y 

110: 

91 b "MUM" 

111: 

X/le3+X 

112: 

Y+r  12 

113: 

if  X # Y ? fto 

" I e r r " 

114: 

fib  "DATA" 

115: 

D+r8 

116: 

f sb  "DATA” 

117: 

D+r9 

118: 

if  r 8 » r 9 5 

eto 

" i nt  food " 

119: 

" I e r r " J 

120: 

fxd  2 

121: 

pet  " T i n i n f 

error!  " 

122: 

prt  "clock 

ticks  = " f r S 

123: 

prt  " inter r 

UPt  s 

*="  > r9 

124: 

prt  " I n t e r v 

al  # 

" » r 1 0 

125: 

prt  "start 

t ine 

: " < T [43 

126: 

prt  "End 

t i ne 

: " ? T [ 1 3 

127: 

p r t " # t i n e 

increnents: " * 

T C33 

128: 

prt  "new 

i nt- , 

(nsec)  = " * 

r 1 2*  1 e3 

129: 

p r t " ** 

Note!"Jprt  "Tin 

i ni 

based  on  # 

of  interrupts" 

130: 

" intfood" : 

131 : 

prt  " 

ii 

132: 

prt  " i n t e r v 

al"  > 

r 125 rl2+Y 

133: 

prt  " 

•• 

134:  src 
135:  O + rll 
136:  TC4] M r 9 - 
2) T C2] + . 15  + TU3 
137:  if  Tt5]*05 
I + TC53  J0+TU3 
138:  T C 1 3 -*T  C43 
139:  Y + T[23 
140:  0 + TI33 
141:  r 10+ 1-M- 10 
142:  fto  "CSI" 
143:  "S": 

14  4:  if  r 1 3 * 0 : 

1 + r 1 3?  0 + T [ 1 3 


145:  fib  M DATA" 
146:  if  2>605 
fto  "Serr" 

147:  if  L«05sto 
" 0 p r i n t " 

148:  fts  (TCI]  ) ■* 
Ct  [L»  0 C L 3 3 
149«  fts  (D)+CJ[ 
L » Q CL 3 +43 
150:  QCL3+8+QCL3 


151: 

ii 

fto  " D p r l n t 

152: 

"Serr": 

153: 

prt-  "Data 

Exponent.  “ 

154: 

prt  " >60" 

155: 

prt.  "Serr " 

156: 

StP 

157: 

fto  "CSI" 

158: 

"Dprint  “ : 

159: 

if  r 11=05 

1 + r 1 1 5 p r t " C 

T i ne 

Data" 

160: 

f n t 8 > f 1 . O t 

e9 . 2 

! » f 6 . 3 

161 : 

w r t 1 6 ■ 8 1 C » 

T [ 1 3 

>D 

162: 

fto  "CSI" 

163: 

" * " : 

164: 

T C 1 3 +T  [23 +T 

[13 

165: 

T [33  +1  + T [33 

166: 

fto  “CSI" 

167: 

" 0 " : 

168: 

: 

169: 

"space " : 

170: 

..  p .. . 

171 : 

" 0 " : 

172: 

fto  "CSI" 

173: 

“csif  in": 

174: 

b e e p 5 e n t 

"Store  on  tape? 

” I sJ  ^ 

175«  if  caP(Jf)# 
"Y"lfto  "ns top" 
176:  enp  "Track 
#"*Y 

177:  enp  "File 
# > X 

178:  trk  YJfdf  X 
179:  idf  Z C 1 3 * 

Z [23  i Z C33  * Z [43  » 
ZC53 

180 J if  Z [23  #05 
tnt  "Overwrite 
f ilt??"i J* 
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181 « if  Z [23  #0  5 
if  cap ( Jt ) # “V" : 
fto  "csifin" 
182:  if  Z [43 <50O 
05be*Pl*nt  "Fil 
e too  snall"« 
J*5fto  "csifin" 
183:  ref  X>C* 
184:  p r t "file 
stored" 

18S:  "nstop" : prt 
"Done " 5 spc  3 5 

StP 

+ 2734 
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